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1 .  This  is  the  Pinal  Report  on  the  contract.  It  covers  the  period 

1  July  77  to  11  Oct  79.  The  objective  of  this  research  was  to  study 
the  waves  which  can  propagate  on  a  honeycomb  panel  of  the  type  used 
to  construct  the  shelters  which  house  much  of  the  ground-based  electronic 
equipment  developed  by  the  Electronic  Systems  Division.  By  knowing 
the  waves  that  can  propagate  on  such  panels,  it  will  be  possible  to 
develop  a  system  of  non-destructive  evaluation  (NDE)  which  can,  in  a 
cost-effective  way,  find  defe  tive  panels  and  avoid  their  incorporation 
into  shelters. 

2.  Although  t^is  work  terminated  before  developing  a  complete  theory 
appropriate  to  the  exact  honeycomb  structures  used  in  shelter  panels, 
it  does  provide  adequate  information  for  understanding  the  waves 
appropriate  for  the  NDE  of  shelter  panels  and  would  form  the  basis  of 
an  experimental  study. 

ALAN  J.  BUDREAU 
Project  Engineer 
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I.  INTRODUCTION 


Honeycomb  panels  and  foam-core  panels  consist  of  a  light-weight 
core  bonded  between  face  plates.  This  fabrication  results  in  a  strong 
light-weight  and  insulating  construction  element.  Such  panels  are  used 
in  containers  for  housing  electronic  equipment.  Honeycomb  panels  are 
also  used  in  airplane  construction,  while  foam-core  panels  are  used  for 
insulating  tanks. 

In  the  case  of  the  honeycomb  panel,  the  core  consists  of  an  array 
of  thin-walled,  hollow  cylinders  of  hexagonal  cross-section,  as  shown 
in  Figure  1,  and  resembles  a  honeycomb.  The  core  material  consists 
of  waxed  cardboard  in  the  case  of  equipment  containers,  while  in  air¬ 
plane  construction  it  is  usually  metal.  Face  plates  are  attached  to  both 
ends  of  the  core  using  epoxy  or  other  adhesive  bonding  materials  (one 
of  the  face  plates  is  shown  in  Figure  1). 

Foam-core  panels  consist  of  a  foamed  plastic,  such  as  styrofoam, 
between  face  plates.  The  foam  can  be  generated  by  interacting  the 
chemical  agents  directly  on  the  face  plates,  in  which  case  a  welded  bond 
exists  between  the  foam  and  plate.  Alternatively  a  prepared  foam  core 
may  be  glued  to  the  face  plates. 

During  fabrication  and  use,  defects  can  occur  within  the  panel,  and 
are  therefore  not  detectable  by  visual  inspection.  One  type  of  defect  is 
a  region  of  a  face  plate  that  is  not  bonded  to  the  core.  Another  defect  in 
honeycomb  panels  is  the  presence  of  water  in  some  of  the  cells,  which 
can  degrade  a  cardboard  core  and  weaken  the  epoxy  or  adhesive  bond 
between  core  and  faceplate.  Other  defects  are  buckled  or  fractured 
cell  walls. 
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Currently,  acoustic  non-destructive  evaluation  (NDE)  techniques 
are  used  to  locate  the  defects  cited  above  in  honeycomb  panels.  An 
acoustic  beam  in  a  fluid  is  directed  at  normal  incidence  on  the  panel  and 
received  at  the  opposite  face.  The  beam  is  then  scanned  in  a  raster 
fashion  over  the  area  of  the  panel.  Unbonded  regions  interrupt  the  trans¬ 
mission  while  water  in  cells  increases  the  transmission.  Since  the  size 
of  the  acoustic  beam  is  on  the  order  of  the  cell  diameter,  the  raster 
scan  inspection  procedure  is  very  lengthy  for  a  large  panel.  A  great 
reduction  in  inspection  time  could  be  achieved  by  propagating  the  acous¬ 
tic  waves  parallel  to  the  face  plates,  rather  than  perpendicular  to  them. 

In  this  way,  one  would  inspect  an  entire  strip  along  the  panel  at  one  time. 
Reinspecting  installed  panels  with  only  one  face  exposed  would  require 
using  waves  that  propagate  along  the  panel. 

A.  Honeycomb  Panels 

Low  frequency  studies  of  propagation  along  honeycomb  panels^ 
have  been  carried  out  with  the  aim  of  NDE  inspection  in  mind.  How¬ 
ever,  it  was  recognized  that  successful  NDE  inspection  would  require 
the  use  of  high  frequencies  for  which  the  acoustic  wavelengths  would 
be  on  the  order  of  the  cell  diameter.  Establishment  of  a  sensitive  test 
procedure  requires  a  knowledge  of  the  propagation  characteristics  of 
the  modes  guided  by  honeycomb  panels.  From  such  knowledge,  one  can 
determine  the  frequency  and  mode  polarization  most  useful  for  inspec¬ 
tion.  These  choices  must  be  made  in  relation  to  the  geometric  design 
and  elastic  properties  of  materials  used  in  the  panel. 

In  this  report  we  first  compute  the  properties  of  waves  propagating 
through  a  honeycomb  of  infinite  extent.  Because  the  honeycomb  is  a 
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periodic  structure  in  two  dimensions  (and  uniform  in  the  third  dimen¬ 
sion),  the  waves  propagating  through  it  are  Bloch  waves,  which  are  anal¬ 
ogous  to  plane  waves  in  an  elastic  continuum.  Each  elastic  field  quantity 
in  a  Bloch  wave  is  the  product  of  a  plane  wave's  exponential  factor,  and  a 
second  function  that  gives  the  variation  of  the  field  quantity  within  a  unit 
cell.  In  the  case  of  a  honeycomb,  this  second  function  is  defined  only 
within  the  honeycomb  walls  if  the  space  within  the  cells  is  assumed 
empty.  The  waves  guided  by  an  assembled  panel  will  be  a  combination 
of  the  Bloch  waves,  just  as  the  waves  guided  by  a  plate  are  a  combina¬ 
tion  of  the  plane  waves  of  an  infinite  medium. 

The  calculation  of  Bloch  wave  propagation  characteristics  are  car¬ 
ried  out  over  the  range  from  low  to  high  frequencies.  At  low  frequen¬ 
cies  (below  1,0  kHz  for  a  Shelter  Core  honeycomb)  the  core  is  found 
to  act  as  an  elastic  continuum  of  hexagonal  symmetry,  and  hence  sup¬ 
porting  three  wave  types.  At  high  enough  frequencies  (above  10  kHz 
for  the  Shelter  Core  honeycomb)  ,  the  wavelength  of  all  wave  types 
is  on  the  order  of  the  cell  diameter,  and  hence  the  periodicity  of  the 
honeycomb  strongly  effects  the  propagation. 

Having  computed  the  Bloch  waves  in  an  infinite  honeycomb,  we  next 
compute  the  characteristics  of  the  waves  guided  by  a  structure  consist¬ 
ing  of  one  face  plate  bonded  to  a  semi-infinite  honeycomb.  In  this  case, 
the  waves  guided  by  the  face  plate  are  strongly  perturbed  by  the  coupling 
to  the  core.  The  SH  and  symmetric  Lamb  waves  in  the  face  plate  are 
found  to  radiate  or  leak  energy  in  to  the  honeycomb,  and  hence  are  of 
the  leaky  wave  type.  The  flexural  wave  in  the  face  plate  is  speeded  up  by 
the  honeycomb,  but  remains  a  bound  surface  wave.  These  calculations 
are  carried  out  treating  the  core  as  an  elastic  continuum.  For  a  Shelter 


Core  honeycomb  this  means  frequencies  below  1.0  kHz  for  the  SH  wave 
in  the  face  plate,  and  below  10  kHz  for  the  symmetric  Lamb  and  flexural 
modes  of  the  face  plate. 

The  cross-section  of  a  regular  hexagonal  honeycomb  is  shown  in 
Figure  2.  The  core  consists  of  strips  or  plates  of  thickness  t,  width  w 
and  infinite  length  along  z  (out  of  plane  of  the  paper).  The  strips  are  as¬ 
sumed  to  be  arranged  in  a  regular  hexagonal  pattern.  This  assumption 
idealizes  actual  honeycombs,  whose  fabrication  leads  to  hexagons  having 
unequal  vertex  angles,  and  two  of  whose  side  walls  are  of  double  thick¬ 
ness  as  compared  to  the  other  four. 

A  unit  cell  of  the  periodic  core  is  indicated  by  the  dashed  curve  in 
Figure  2.  The  entire  periodic  core  can  be  generated  by  translations  of 
the  unit  cell  along  the  basis  vectors  d^  and  d^,  as  shown  in  Figure  2, 
of  the  periodic  structure.  The  cell  is  seen  to  be  composed  of  plate 
segments  coupled  at  two  Y  joints.  A  primary  part  of  the  study  of  core 
waves  therefore  involves  the  coupling  of  the  plate  waves  at  the  Y  joint. 

The  acoustic  fields  in  the  plate  segments  that  form  the  cell  walls 
can  be  expressed  as  a  superposition  of  plate  modes.  Since  the  thickness 
t  of  the  cell  walls  is  one  tenth  of  their  width  w,  or  less,  thin  plate  ap¬ 
proximations  may  be  made  for  the  plate  modes,  even  for  frequencies 
where  the  cell  diameter  is  on  the  order  of  the  acoustic  wavelength. 
Consistent  with  this  approximation  we  keep  only  the  lowest  four  plate 
modes,  which  are:  1)  the  shear  wave  with  no  variation  across  the  plate  - 
the  SH  mode;  2)  the  lowest  symmetric  Lamb  mode  -  L  mode;  3)  the  lowest 
anti-symmetric  Lamb  mode  or  flexural  wave  -  F(0)  mode;  and  4)  the  first 
cutoff  flexural  wave  -  F(l)  mode. 
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gure  2.  Cross-section 


The  F(l)  mode  is  included  since  it  has  the  same  transverse  varia¬ 
tion  through  the  thickness  of  the  plate  as  the  F(0)  mode,  and  since  its 
attenuation  constant  vanishes  as  frequency  goes  to  zero.  Hence,  this 
mode  will  couple  neighboring  Y  joints  at  low  frequency  and  its  inclusion 
is  necessary  to  accurately  find  the  amplitude  of  the  F(0)  mode  excited  even 
at  an  isolated  Y  joint.  All  higher  order  plate  modes  attentuate  rapidly 
away  from  the  Y  joint,  even  in  the  limit  of  zero  frequency,  and  there¬ 
fore  do  not  couple  neighboring  Y  joints.  In  addition  these  modes  have 
more  variation  through  the  thickness  of  the  plate,  and  are  therefore  not 
expected  to  be  strongly  excited. 

The  physical  approximations  used  to  obtain  the  mode  coupling  are 
discussed  in  Section  II.  These  approximations  are  embodied  in  a  set  of 
twelve  equations  containing  the  amplitudes  of  the  twelve  modes,  (four  modes 
in  each  of  the  three  plates).  From  these  equations  a  scattering  matrix  is 
obtained  that  gives  the  amplitudes  of  each  of  the  twelve  modes  excited  at 
the  Y  joint  for  any  incident  mode. 

For  propagation  perpentidular  to  the  cell  axes  the  SH  mode  does  not 
couple  to  the  L,  F(0)  and  F(l)  modes  at  the  Y  joint,  as  discussed  in  Sec¬ 
tion  II.  Thus  we  may  treat  separately  the  case  of  waves  having  particle 
motion  parallel  to  the  cell  axis,  which  corresponds  to  the  SH  plate  mode 
polarization,  and  whose  properties  are  found  in  Section  III.  The  L,  F(0) 
and  F(l)  modes  together  give  rise  to  two  Bloch  waves  in  the  honeycomb 
polarized  with  particle  motion  perpendicular  to  the  ceil  axis,  which  are 
considered  in  Section  IV. 

The  case  of  propagation  oblique  to  the  cell  axes  is  treated  in  Section 
V.  Here,  the  shape  of  the  intersection  of  the  slowness  surfaces  with  the 
principle  planes  is  found.  From  these  curves  and  the  propagation 
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perpendicular  to  the  cell  axes  it  is  possible  to  infer  the  shape  of  the  com¬ 
plete  slowness  surface.  At  high  frequencies  it  is  found  that  more  than 
three  surfaces  can  exist,  and  these  surfaces  indicate  that  propagation 
can  take  place  only  over  a  limited  range  of  angles  about  the  direction  of 
the  cell  axes. 

The  characteristics  of  Bloch  wave  propagation  at  low  frequencies 
are  used  to  define  the  constants  of  an  elastic  continuum  of  hexagonal  sym¬ 
metry.  This  continuum  is  assumed  to  be  directly  bonded  to  an  aluminum 
face  plate  and  the  propagation  constant  of  waves  guided  by  the  combina¬ 
tion  are  determined  in  Section  VI. 

B.  Foam-Core  Panels 

In  studying  foam-core  panels,  we  have  treated  the  case  of  waves 
guided  by  a  single  aluminum  face  plate  attached  to  a  semi-infinite  foam 
core,  as  discussed  in  Section  VII.  The  foam  is  treated  as  an  isotropic 
continuum,  with  elastic  constants  representative  of  typical  static  values 
found  in  the  literature.  Both  the  case  of  directly  bonded  and  glued  foams 
are  considered.  The  waves  are  found  to  radiate  or  leak  energy  into  the 
foam.  The  rate  of  leakage  and  the  propagation  constant  are  found  to  be 
strongly  influenced  by  the  presence  of  a  glue  layer. 


II.  MODE  COUPLING  AT  A  Y  JOINT 


The  Y  joint  of  three  plates  is  shown  in  Figure  3.  For  each  plate, 
with  index  i  =  1,  2,  3  a  coordinate  system  (x^,  y^,  z)  is  established  as 
shown.  In  treating  Block  wave  propagation  oblique  to  the  cell  axes, 
all  four  of  the  plate  modes  in  the  cell  walls  will  couple  at  the  Y  joint. 

In  order  to  satisfy  the  coupling  equations  for  all  values  of  z  along  the 
joint,  all  of  the  modes  must  have  the  same  phase  dependence  expj-jk^z), 
where  k  is  the  wavenumber  along  z.  This  condition  is  the  analogue 
of  Snell's  law  in  optics.  Since  the  Bloch  wave  is  formed  by  the  plate 
modes  that  are  multiply  scattered  at  the  Y  joint,  all  fields  in  the  Bloch 
wave  must  have  the  same  phase  variation  along  z,  and  hence  k  is  the 
z  component  of  the  wavevector  of  the  Block  wave.  The  variation  of  the 
modal  fields  along  (x.,  y.,  z)  in  each  plate  can  be  represented  by  equi- 
valent  transmission  lines'  .  Mode  coupling  at  a  Y  joint  is  then  repre¬ 
sented  by  an  interconnection  of  the  modal  transmission  lines,  or  by  a 
scattering  matrix. 


A.  Transmission  Line  Representation 

In  the  transmission  line  representation,  the  fields  of  each  mode  for 
harmonic  time  dependence  exp(jwt)  are  expressed  in  the  form 
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In  (1),  the  1x3  matrices  £(y^)  and  cjjy^)  are  the  mode  vectors  and  V(x^), 
I(x^)  are  the  modal  voltages  and  currents.  In  terms  of  traveling  waves. 
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where  K  is  the  wavenumber  of  the  mode,  Z  is  the  modal  impedance,  and 
b,  a  are  the  amplitudes  of  the  waves  traveling  in  the  +  and  -  x.  directions, 
respectively.  If  the  mode  vectors  are  normalized  such  that 
t/2 
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then  Re[V(x.)I  (x.)j  represents  the  power  carried  in  the  +  x.  direction 
by  the  mode,  per  unit  length  along  z.  For  propagating  modes  K  and  Z 
are  real  and 


Re[V(x.)I*{x.)J  =  |  b|  2  -  ja|2. 
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while  for  evanescent  modes  with  K  imaginary,  Z  is  imaginary  and 

Re[V(x.)I*(x.)]  =  Re^y^j-  (ab*  -  a*b)j  .  (5) 

The  various  modal  quantities  for  plate  waves  have  previously  been  de¬ 
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the  plane -strain  approximations  may  be  made  for 
In  this  approximation,  the  cross-section  of  the 


plate  is  assumed  to  remaiu  a  plane  surface  of  constant  width.  Employ¬ 
ing  this  approximation,  we  derive  in  the  Appendix  the  various  quantities 
for  each  mode.  The  modal  quantities  are  tabulated  below,  in  these 
pression.,  u  is  the  shear  modulus  of  the  wall  material,  p  is  it,  den- 
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and  r  is  given  by  (13), 

B,  Representation  of  Coupling  at  Y  Joints 

The  Y  joint  of  three  plates,  as  shown  in  Figure  3,  may  be  viewed 
as  the  connection  to  a  triangular  cylinder  of  three  plates  having  end 
faces  perpendicular  to  x^.  With  this  view,  mode  coupling  at  the  Y 
joint  is  determined  by  the  boundary  conditions  at  the  surfaces  of  the 
triangular  cylinder.  We  assume  that  the  cylinder  is  rigid  and  mass¬ 
less.  The  assumption  of  rigid  motion  composed  of  translation  and 
rotation  about  the  z  axis  is  consistent  with  the  plane  strain  approxi¬ 
mation  used  for  the  modal  fields.  Neglecting  the  mass  of  the  cylin¬ 
der  implies  that  the  net  force  and  torque  about  the  z  axis  must  vanish. 


and  ignores  the  kinetic  energy  that  it  stored  in  the  cylinder's  motion. 
Rotation  of  the  cross-section  about  the  x  and  y  axes  is  not  taken  into 
account,  nor  are  the  torques  about  these  axes  considered.  Experience 
with  structures  that  are  periodic  in  one  dimension  suggests  that  the  ef¬ 
fect  of  neglecting  the  kinetic  and  potential  energy  in  the  cylinder  is 

limited  to  a  small  error  in  the  frequencies  predicted  for  the  stop- 

(9) 

bands  ,  which  is  not  expected  to  be  significant  for  NDE  applications. 
The  approximations  used  here  are  essentially  the  same  as  those  used 
to  describe  mode  coupling  at  E  and  butt  joints  of  two  plates^ 

If  the  rotation  and  torque  about  the  x  and  y  axes  were  taken  into 
account  more  equations  would  be  obtained  than  the  number  of  scat¬ 
tered  waves.  As  a  result,  the  scattering  problem  would  be  over 
determined,  and  it  would  be  necessary  to  include  more  scattered 
modes  in  the  plates.  Physically,  the  particle  velocities  associated 
with  the  rotation  of  the  cross-section  about  x  and  y  will  be  small 
compared  to  the  translational  velocities.  To  estimate  the  differences 
in  the  velocities  cited,  consider  the  F(0)  mode  to  be  incident  in  plate 
1  of  Figure  3  at  an  angle  0  with  respect  to  the  x^  axis.  This  mode 
will  tend  to  rotate  the  triangle  about  the  axis,  in  addition  to  rotat¬ 
ing  it  about  z  and  translating  it  along  y^.  At  the  corner  x^  =  0, 
y^  =  t/2,  the  particle  velocity  due  to  the  rotation  about  x^  is  v^(t/2) 
while  the  translational  velocity  is-  v^.  From  the  mode  vector  given 
by  (21)  it  can  be  shown  that  the  ratio  vz(t /2)/v^  is  of  the  order 

/  k  t/2  sin  9,  which  is  less  than  about  (0.  i)sin  6  for  the  highest 
'a 

frequencies  discussed  here.  In  other  words,  the  particle  velocity 
due  to  rotation  about  x^  is  small  compared  to  the  translational  velocity. 
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and  hence  little  error  will  be  introduced  by  not  accounting  for  the 
rotation. 

1 .  Particle  Velocity  Conditions 

Let  v^(y^)exp(-jkzz)  be  the  particle  velocity  at  the  end  face  x^  =  0 
of  the  ith  plate  in  Figure  3.  Now  the  motion  of  the  triangular  cylinder 
consists  of  a  translation  v  exp(-jk  z)  of  the  center  of  gravity  plus  a 
rotation  ft  z_0exp(-jk  z).  Referring  to  Figure  4,  it  is  seen  that  the 
particle  velocity  of  the  ith  face,  with  the  phase  factor  exp(-jkzz)  un¬ 
derstood,  is 
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where  x  .  and  y  .  are  unit  vectors  along  x.  and  y.,  and 
— oi  -*-'01  B  l  7 1 

t 


d  = 


2  fT 


(19) 


is  the  perpendicular  distance  from  the  center  of  gravity  to  the  trian¬ 
gular  faces. 

Alternatively,  the  particle  velocity  Vj(y  )  can  be  expressed  in 

terms  of  the  modes  of  the  plate.  Let  Vgj,  Ig^  be  the  voltage  and 

current  at  x.  =  0  of  the  SH  mode  in  the  ith  plate.  Similarly,  let 

VT  .,  IT  .  -and  Vn.,  I„.  and  V.  T  .  be  the  voltages  and  currents  at 
Li  Li  Oi  Oi  li  li  6 

x.  =  0  of  the  L,  F(0)  and  F(l)  modes  in  the  ith  plate.  The  particle 
velocity  at  x.  =  0  can  be  expressed  in  terms  of  the  plate  modes  using 
the  modal  quantities  listed  in  the  previous  section.  With  the  phase 
factor  exp(-jk  z)  implied,  the  components  of  particle  velocity  are 
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Figure  4 


Cross-section  of  the  triangular  cylinder  at  a  Y  joint,  which 
is  assumed  to  move  as  a  rigid  body. 
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In  writing  (22)  the  z  components  of  particle  velocity  of  the  F(0)  and 
F(l)  modes  were  not  included  since  they  depend  linearly  on  y^,  and 
hence  represent  a  rotation  about  the  x^  axis,  which  is  not  accounted 
for  in  our  approximation. 

The  particle  velocity  (20)-(22)  of  each  plate  end  must  be  consistent 
with  translation  of  the  triangular  cylinder  and  its  rotation  about  z,  as 
stated  in  equation  ( 1  8).  Equating  like  powers  of  y^  for  each  component 
gives  the  equations 
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Four  equations  can  be  obtained  by  substituting  (23)  into  (25)  and 
using  the  relations 
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-O'  2  Xol  ‘  v 


Z°2‘  -°=  "  ^fzol  ■  Zo-^Lol  •  v, 


(27.  c) 


Xoi‘2o-  1TL £oi-£«-|i01- 

A.so,  repress  two  e,uations  reUting  the  modsi 
<ioeM26,.  Th.  eight  e,ua,io„a  so  „btlined  ar„ 

°  ”  7W  <Il2  +  2  'w  ’  ■  j=-lT  <!S2  +  |  ISl ) 


(27.  d) 


[No^-U-N7^j+jM|i5^+fJ_U_  (23.a) 


0  *  /IT  +  2  lU  ' 

-  f  V°1  Vu  ]  f  I„  I 

fir*  LNo,2r-‘  ‘  ■  iksd  foi^n +  n-tM 


^21+1)  <28.  b) 


(28.  c) 


0=  ^=r  It.!  If-L-floiliZoi  V13+^V 


+  jk_d  &lgjoi  iI13+iIiil 

8  p^r-i  r +  i^rs+rr  j 


(28.  d) 


-19- 


*01  ,  *11  l02  *12 
NQ( 2r-i )  +  N1(2r  +  i)  "  NQ(2r-l)  N1(2r  +  1) 

-  X03  X13 

"  N0(2r-1)  +  N1(2r+l) 


(1  ~V)  "TX  VLi  +  VS1  =  (i  ~V)  ~T  VL2  +  VS2 


=  VL3  +  VS3 

ZkL 


2.  Force  and  Torque  Conditions 

In  addition  to  imposing  continuity  of  velocity,  which  yielded  the 
foregoing  six  equations,  we  must  require  that  the  force  and  torque 
about  z  on  the  cylinder  vanish.  A  mathematical  statement  of  vanish¬ 
ing  force  is 
3  t/2 

£  [^1  2  Txx^0,  yi*  z)  +  ^°i  Txy(0*  yi*  z)  +  — °  Txz(0,  yi’  dyi  =  0 


while  the  vanishing  of  torque  about  the  z  axis  can  be  written 
3  t/2 

l  f  (*'old  +  lol  y i>  x  [iol  Txx(0'  yi’  z)  +  Zol  Txy(0'  yi’  z)]  dyi  =  ° 
i=l  -t/2 


Substituting  the  modal  quantities  from  the  previous  section  into 


(31)  gives 
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k  \  k 

t+d-^) -§  \ft »L.f!iS  /Tv 
kLy  8 


^■V1  Oi  ‘,i 


+  *oI  1  2(1  -1/  )  l  N0  ■  Nj 


k  2k  ^  -  k2  "1 

+  2„  /mm  -2iTILi+-^T-^  [Si  }  =  ° 
-  Ks  -1  J 


Decomposing  x  y  for  i  =  2,  3  along  x  ,  y  we  can  write  (33) 
— Oi  oi  — oi  oi 

in  component  form  as 


k  k 

o  =  -l+«  -V)  -§  ft  (Vu  -  i  vL2  -  i  vL3)  +  /T  (2VS1-VS2-VS3) 

klJ  s 


Si  S2  S3 


^'kst  /^102~I03  I12~113 

tfT(i-v)  \  No  '  Ni  , 


^  (VI,J-VL3 


<V,  -  V,>  +  /3*  r5-  (v<=,  -  v.,1 


k  v  S2  S3' 
s 


/^~ks  *  /2I01  ~  *0  2  ~  *03  2I1 1  “  1i  2  "  *1  3 


24(l-i/) 


k  2k2-k2 

0  =  “2  k”  (IL1  +  JL2+  IL3)  +  (ISi  +  *S2  +  *83* 

s 

Substituting  the  modal  quantities  into  (32)  gives 
°*  £ ks<(^ ' +  ^  ( '2 ^ 
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3.  Scattering  Matrix  Representation 


Equations  (28)-(30)  and  (34 ) -(3 7 )  are  twelve  equations  describing 
the  coupling  between  the  twelve  plate  modes,  four  in  each  plate.  Be¬ 
cause  of  their  complexity,  we  have  chosen  to  use  the  equations  to 
generate  a  scattering  matrix  for  the  junction.  To  this  end  we  express 
the  various  voltages  and  currents  in  terms  of  the  traveling -wave  am¬ 
plitudes  a  and  b  using  (2)  with  x^  =  0. 

Let  the  1x4  matrices  a.  and  b.  be  defined  as 


aSi 

bSi 

aLi 

kL  = 

bLi 

a0i 

b0i 

ai  i_ 

.v 

(38) 


After  substantial  manipulation,  one  obtains  from  (28)-(30)  and  (34)- 
(37)  equations  relating  the  traveling -wave  amplitudes  in  the  form 


(39) 


The  12x12  matrices 

r  (-)l 

and  A;.  are  defined  in  Table  I.  From 

L  ij  j 

L  ij  J 

(39).  the  1  2x1  2  scattering  matrix  relating  the  a's  and  b's  as 


D>]  =  [s,j]  W 


(40) 


can  be  expressed  in  the  form 
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TABLE  I 


Conditions  (42)  for  the  modal  voltages  and  currents  can  be  interpreted 
readily  in  terms  of  the  interconnection  of  the  equivalent  modal  trans¬ 
mission  lines.  Since  the  voltage  is  the  same  on  all  lines  and  the  cur¬ 
rents  sum  to  zero,  this  transmission  line  connection  is  as  shown  in 
Figure  5,  and  is  used  subsequently  to  find  the  Bloch  waves  propagat¬ 
ing  in  the  honeycomb  and  polarized  with  particle  motion  along  z. 

While  the  coupling  of  the  L,  F(0)  and  F(l)  modes  is  also  simpler 
for  k  a  0,  it  is  still  not  obvious  or  useful  to  introduce  equivalent  cir- 
cults  interconnecting  the  equivalent  transmission  lines.  In  this  case 
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it  is  still  convenient  to  use  the  scattering  matrix  formation,  where 
and  are  now  1x3  matrices  obtained  by  dropping  ag^  and  bg.  in 
(39).  Also  j  and  ^A^j  are  9x9  matrices,  whose  elements  in 

Table  II  can  be  obtained  by  dropping  the  first,  fifth  and  ninth  columns, 
and  seventh,  eighth  and  eleventh  rows  of  Table  I.  Also,  the  order  of 
the  rows  has  been  changed,  which  does  not  effect  the  value  of  the 
scattering  matrix. 
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III.  BLOCH  WAVES  POLARIZED  ALONG  THE  CELL  AXES  FOR  PROPAGATION  PERPEN¬ 
DICULAR  TO  THE  CELL  AXES 


Because  the  SH  plate  waves  do  not  couple  to  other  modes  at  the  Y  junc¬ 
tions,  and  because  they  are  non-dispersive ,  the  Bloch  waves  composed  of 
them  have  the  simplest  properties.  We  therefore  consider  this  polariza¬ 
tion  first,  in  order  to  demonstrate  the  effects  of  the  six-fold  symmetry 
on  propagation  characteristics. 

A.  Formal  Solution  for  Dispersion  Characteristics 

Starting  with  the  equivalent  circuit  for  the  Y  junction,  as  depicted 
in  Figure  5,  the  equivalent  network  for  a  unit  cell  can  be  obtained.  This 
network  is  drawn  in  Figure  6  and  is  seen  to  consist  of  four  ports  at  which 
the  voltages  and  currents  ,  1^  (i  3  1  -  4)  are  defined.  The  voltages  and 
currents  can  be  related  via  an  impedance  matrix.  Because  of  the  symmetry 
of  the  network,  this  relationship  takes  the  form 


"V 

V2 

V„ 

3 

V4 

— 

11 

Z12 

Z13 

Z13 

12 

zu 

Z13 

Z13 

13 

Z13 

Z11 

Z12 

13 

Z13 

Z12 

Z11 

(43) 


By  straightforward,  but  lengthy  network  analysis,  it  can  be  shown  that 
Zn’  zi2 and  zi3  are  siven 


7  .  47  1-10  tan28  +  7  tan^Q 

^11  ^  s  2 

6  tan  9  (2tan  0-1) 


(44) 


—  jZ  o  9 

6  cos40  tan  0  (2tan40-l) 


1-5  tan  0 


(45) 


'13 


JZ. 


6  cos40  tan  9  (2tan  0-1) 


(46) 


with 


0  *»  k  w/2 
s 


(47) 


We  now  consider  the  terminal  voltages  and  currents  and  1^  to  be 


associated  with  a  wave  traveling  through  the  core.  Let 

k  =  x  k  +  y  k 
—  —ox  4)  y 


(48) 


be  the  wave  vector  of  this  wave.  Also,  let  d^  and  d^  be  the  basic  trans¬ 


lation  vectors  of  the  periodic  structure,  as  shown  in  Figure  2,  and  given  by 

(49) 


d,  =  p  w  x  +  w  y. 
-1,2  2  -o  2  ^o 


The  Floquet  assumption  for  periodic  structures  then  implies  the  following 
relation  between  the  terminal  voltages  and  currents 

-jk  •  d. 


V3  =  Vx  e 


-1 


-jk 


i3  =  -  Tl  e 


-1 


(50) 


V4  ~  V2  6 


-jk  *  —2 


h  6 


-jk  •  d2 


The  minus  sign  before  1^  and  I2  is  due  to  the  assumption  that  the  currents 
are  positive  into  the  unit  cell. 

Substituting  (50)  into  (43)  one  obtains  four  linear  homogeneous  equa¬ 
tions  in  terms  of  the  four  unknowns  V^,  V^,  1^  and  1^.  Defining 


^  “  e 


■  jk  * 


-jk  •  d2 


(51) 
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the  four  equations  may  be  written  as 


Z11  '  ^ 1Z13  Z12  ‘  *2Z13  _1  0 

*1 

0  = 

Z12  ‘  ^1Z13  Z11  '  ^2Z13  °  "1 

T2 

Z13  "  ^  1Z1 1  Z13  ‘  ^ 2Z 12  -*1  0 

V1 

Z13  *  *1Z12  Z13  "  *2Z11  °  "*2 

V2 

— 

L  J 

(52) 


In  order  for  (52)  to  have  a  non-trivial  solution,  the  determinant  of 
the  4x4  matrix  must  vanish.  Setting  the  determinant  equal  to  zero  gives 
the  dispersion  equation  of  the  form 


(Z13_Z12)COS  6  +  2Z1,(Z1,-Z11)cos  6  cos  y  +  (Z^-Z,,)  =  0 


"13  12  "11 


11  13' 


(53) 


where 


(54) 


This  equation  may  then  be  solved  numerically  for  its  root  k^  as  a  func¬ 
tion  of  k  and  m  or  k  as  a  function  of  k  and  it.  The  root  corresponds 

y  y  * 

to  a  wave  in  the  core  whose  polarization  is  found  from  the  eigen  vector 
associated  with  the  root. 

Because  of  the  periodicity  of  the  honeycomb,  (53)  will  have  multiple 
solutions.  Let  u^  and  u^  be  reciprocal  lattice  vectors 


~1>2 


=  2n  x  +  2n _ y 

3w  -o  +  /3  w  *o 


(55) 


so  that 


^ij 


(56) 
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where  6_  is  the  Kronecker  delta.  Thus,  if  a  particular  wavevector  k 
satisfies  (53),  then  for  any  integers  m  and  n  the  wavevector  k  +  m  u^  + 
n  u^  will  also  satisfy  (53).  Also,  due  to  the  hexagonal  symmetry  of  the 
honeycomb,  a  plot  of  k^  versus  k^  for  fixed  uu  will  have  six-fold  symmetry 
about  the  origin. 

B.  Numerical  Evaluation  of  Propagation  Characteristics 

Numerical  computations  have  been  carried  out  for  the  dispersion  curves 

k  versus  k  for  various  fixed  values  of  u).  Because  the  SH  wave  velocity 
x  y 

is  independent  of  the  wall  thickness  t,  and  because  the  impedances  of  all 
transmission  lines  in  Figure  6  are  the  same,  the  results  are  independent 
of  t.  In  addition,  if  the  normalized  frequency  variable  8  in  (47)  is  used, 
and  the  wave  numbers  k^  and  k^  are  multiplied  by  w/2,  the  computations 
will  apply  to  all  materials  and  honeycombs. 

Due  to  the  symmetry  and  periodicity  of  the  honeycomb,  the  first  quad¬ 
rant  in  the  k  -k  plane  of  the  first  Brillouin  zone  is  more  than  adequate 
x  y  M 

to  present  the  numerical  results  for  the  dispersion  curves  of  k^  versus 
ky.  However,  in  order  to  explain  the  variation  of  the  dispersion  curves 
with  frequency  we  have  sketched  the  curves  for  several  Brillouin  zones  in 
Figures  7(a),  (b),  (c).  In  Figure  7(a)  we  have  indicated  several  reciprocal 
lattice  points  by  dots  and  the  vectors  u^  and  u^.  The  dashed  hexagon 
represents  the  first  Brillouin  zone.  For  low  frequencies  (8  «  0.6155), 
the  dispersion  curve  is  very  nearly  a  circle  of  radius  |k  |  ^  ©s. 

This  value  of  |k|  corresponds  to  a  wave  velocity  for  the  Bloch  wave  that 
is  l/ypl  times  the  shear  wave  velocity  of  the  honeycomb  material.  Because 
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i 


ro|£ 


of  the  periodicity,  the  near-circles  centered  about  each  reciprocal 
lattice  point  are  all  locii  of  solutions  of  the  dispersion  relation,  (53), 
as  discussed  after  (56). 

As  the  frequency  variable  9  increases  to  just  below  0.6155,  the  near¬ 
circular  dispersion  curves  expand  and  distort  into  nearly  hexagonal  shapes, 
as  indicated  in  Figure  7(b).  At  9  =  0.6155,  the  tips  of  the  hexagons 
touch.  For  0.6155  <  9  <  ^  =  0.7854,  the  dispersion  curves  switch  to  being 
nearly  triangular  and  centered  about  the  apexes  of  the  dashed  hexagon,  as 
shown  in  Figure  7(c).  As  9  approaches  rr/4,  the  dispersion  curves  shrink 
to  the  points  at  the  apexes  of  the  hexagon.  For  9  increasing  above  tt/4, 
the  change  in  the  dispersion  curves  is  reversed  from  that  described  above 
until  for  0  approaching  rr/2,  the  curves  become  nearly  circular  as  in  Figure 
7(a).  For  0  increasing  above  rr/2,  the  variation  of  the  dispersion  curves 
described  above  is  repeated  periodically  with  period  tt/2. 

The  actual  dispersion  curves  in  the  first  quadrant  of  the  first  Brillouin 
zone  are  plotted  to  scale  in  Figure  8,  and  are  seen  to  be  in  agreement  with 
the  discussion  of  Figure  7.  One  of  the  features  of  these  curves  that  may 
be  important  for  NDE  is  the  fact  that  in  certain  frequency  ranges  the  waves 
are  cut-off  along  +  x  and  the  four  other  directions  making  angles  of  60° 
to  the  +  x  axes.  These  cut-off  directions  are  ones  parallel  to  a  set  of 
cell  walls.  However,  the  waves  are  never  cut-off  along  y,  and  the  six¬ 
fold  related  axes. 

In  order  to  show  the  cut-off  behavior  for  propagation  along  x,  we  have 
plotted  as  a  function  of  id  for  the  case  k^  ■  0.  This  plot  is  shown  in 
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Figure  8.  Slowness  curves  in  the  first  quadrant  of  the  first  BriHouin 
zone  of  the  SH  Bloch  wave  for  several  values  of  0  =  k  w/2. 


Figure  9,  where  is  real  for  0  <  Q  <  0.6155  and  0.9553  >  6  >  rr/2. 

In  the  range  0.6155  <  0  <  0.9553,  k^  is  complex  with  real  part  given  by 
w 

Re  Gcx  “11/3.  The  curves  are  symmetric  about  0  =  tt/4.  For  9  outside 

the  range  plotted,  the  curves  repeat  with  period  tt/2.  By  contrast,  the 

w  w 

plot  in  Figure  10  of  k^  —  versus  the  frequency  variable  k^  —  for  propaga¬ 
tion  along  y  does  not  show  the  presence  of  a  stop  band.  The  forward  and 
backward  wave  regions  connect  directly  at  6  =  n/4.  Again,  the  curves  as 
a  function  of  0  repeat  with  period  of  n/2. 


The  variations  in  the  magnitudes  of  particle  velocity  and  stress  within 
the  unit  cell  is  depicted  in  Figure  11  for  propagation  along  x  (k^  =*  0) 
and  at  a  frequency  just  below  cut-off.  Because  the  fields  are  symmetric 
about  the  center  of  the  unit  cell,  it  is  only  necessary  to  plot  the  fields 
over  one  half  of  the  cell.  From  Figure  11  it  is  seen  that  the  center  of 
the  unit  cell  is  almost  stationary  but  has  a  large  stress,  as  compared  to 
points  at  the  edges  of  the  cell,  for  propagation  along  x.  At  the  joint, 
however,  the  stress  on  the  center  branch  is  close  to  the  maximum  value  in 
the  cell.  Note  that  the  jump  in  stress  at  the  joint  is  due  to  the  fact 
that  stress  on  the  center  branch  is  the  sum  of  the  stress  on  the  two 
joining  branches,  and  hence  double  that  on  Branch  1  for  propagation  along 
x.  The  high  stress  at  the  joint  would  facilitate  NDE  inspection  for 
cracking  at  the  joint,  since  a  cracked  joint  would  require  stress  to  vanish 
and  hence  would  be  a  large  perturbation.  For  frequencies  well  below  cut¬ 
off,  Iv^  |  and  |tzx  |  are  nearly  uniform  over  the  unit  cell. 
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IV.  BLOCH  WAVES  POLARIZED  PERPENDICULAR  TO  THE  CELL  AXES  FOR  PROPAGATION 
PERPENDICULAR  TO  THE  CELL  AXES 

As  discussed  in  Section  II,  waves  propagating  the  plane  perpendicular 
to  the  axis  and  polarized  with  particle  velocity  in  the  plane  are  composed 
of  L(0),  F(0)  and  F(l)  modes  in  each  cell  wall.  Because  of  the  number  of 
modes  and  the  complexity  of  the  equations  describing  their  coupling  at  a 
Y  joint,  it  is  not  readily  apparent  how  to  represent  the  coupling  by  an 
equivalent  circuit,  as  was  done  for  the  case  of  particle  motion  along 
the  cell  axis.  Instead,  the  scattering  matrix  representation  (40)  is  used 
for  the  coupling.  This  representation  is  convenient  since  the  necessary 
matrix  manipulations  can  be  indicated  symbolically,  whereas  the  actual 
evaluation  is  carried  out  numerically  by  a  computer. 


A.  Formal  Solution  for  the  Propagation  Characteristics 
In  order  to  facilitate  the  matrix  manipulations,  define  the  1x3  column 
matrices 


where  the  a's  and  b's  are  the  traveling  wave  amplitudes  in  the  i  =  1,  2,  3 
plate  of  Figure  3  evaluated  at  x ^  *  0.  Recall  that  the  subscripts  L,  0,  1 
refer  to  the  L,  F(0)  and  F(l)  modes  and  that  the  a's  are  the  amplitudes  of 
waves  traveling  in  the  -x^  direction  towards  the  junction,  while  the  b's 
are  the  amplitudes  of  waves  traveling  in  the  opposite  direction.  With 
this  definition,  the  scattering  matrix  S  may  be  partitioned  into  3x3 
matrices  S.  as 


The  three-fold  symmetry  of  the  Y  junction  is  evident  in  the  form  (58)  of 
S.  In  addition,  all  are  symmetric  and  the  two  matrices  and  are 
the  same  except  for  the  signs  of  the  12,  13,  21  and  23  elements-  Because 
the  F ( 1 )  mode  is  cut-off  rather  than  propagating,  S  is  not  a  unitary  matrix. 
1.  Scattering  Matrix  for  a  Unit  Cell 

Figure  12  shows  a  unit  cell  divided  into  two  identical  portions,  each 
consisting  of  three  plates  of  length  w/2  connected  at  a  Y  joint.  The 


numbering  of  the  plates  has  been  chosen  for  convenience  in  subsequent 
analysis.  If  a^ '  and  b^ '  are  the  traveling  wave  amplitudes  for  the  ith 
plate  at  the  end  x^  =  w/2,  then 


(59) 


where  e  is  the  3x3  diagonal  matrix 


(60) 


(61) 
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from  (59)  it  follows  for  the  left-hand  half  of  the  unit  cell  that 


“  - 

bi 

—2 

=  S' 

—2 

*5 

—5 

(62) 


Similar ly,  for  the  right-hand  half  of  the  unit  cell 


- 

-  — 

b' 

a’ 

-3 

-3 

=  s' 

*6 

(63) 


The  two  halves  of  the  unit  cell  are  coupled  at  the  mid-plane  T,  where 
the  boundary  conditions 


X5X5 


=  T 


V« 


X5y5 


=  T 


X6y6 


(64-a) 


(64-b) 


must  hold.  From  (1),  (2)  and  (10),  the  first  condition  in  (64-a)  and  the 
first  condition  in  (64-b)  are  seen  to  imply  for  the  L  mode  that 


<bLs  +  aLs>  =  <bi,6 +  aLe> 


(b|5 


a^5>  =  '  <bi.6  -  a^- 


(65) 


since  and  q^  are  the  same  for  both  plates.  Because  the  wave  propagating 

in  the  +■  x,  direction  on  plate  5  becomes  the  wave  propagating  in  the  -  x, 
o  6 

direction  on  plate  6,  and  vice  versa,  (65)  is  satisfied  if 


-43- 


'L5  L6 
>’  =  b’ 


(66) 


L5 


L6  J 


When  considering  the  flexural  mode  F(0),  the  second  conditions  in  (64-a)  and 
(64-b)  are  seen  from  (1),  (2)  and  (12)  to  give 


(b05  ‘  a05*  =  (b06  "  a06) 

+  a05>  "  -  <b06  +  a06>- 


(67) 


since  and  are  the  same  for  both  plates.  However,  the  first  conditions 
in  (64-a)  and  (64-b)  are  seen  from  (1),  (2)  and  (12)  to  give 


^05  +  a05)y5  ~  (b06  +  a06)y6 
(b05  '  a05)y5  =_(b06  '  a06)y6 


(68) 


Noting  that  y,-  =  -  y^,  it  is  seen  that  (67)  and  (68)  imply  that 


05 


*06  n 


(69) 


*05  '  b06  J 


where  it  is  recalled  that  the  +  x^  traveling  wave  on  one  plate  becomes  the 
-  x^  traveling  wave  on  the  other  plate. 

Similarly,  for  the  F(l)  mode 


bi5  "  -  ai6 


Jis 


(70) 


'16- 


Equations  (66),  (69)  and  (70)  may  be  summarized  as 


—5 


—  6 


T*6 


T  £5 


(71) 
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where 


T  = 


1 

0 

0 


0 

-1 

0 


0 

0 

-1 


(72) 


In  order  to  join  the  scattering  representations  (62)  and  (63)  using  the 
conditions  (66),  (69)  and  (70)  to  obtain  the  scattering  matrix  for  the  unit 
cell,  we  write  (62)  and  (63)  in  the  form 
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~2 

—3 
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£3 

H 

0 

0 

£3 

si 

i. 

*4 

0 

~2 

-  - 

- 

— 

—  - 

(73) 


and 


(74) 


Substituting  (71)  into  (74),  we  can  solve  for  a £  and  a£  in  terns  of  a|, 
£ 2 ’  £3  and  £4'  If  this  result  is  substituted  into  (73)  we  obtain  the 
scattering  representation  for  the  unit  cell  as 


(75) 
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where  u  is  the  12x12  matrix  given  by 


u  = 


Si  S' 

si  si 

~3  ~1 


0 

0 


0 

0 


0  0 
0  0 


Si  -§2 

^3  ^i 
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s ' 
~2 

0 

0 


*3 

*2 


- 

-1 

—  — 

_si  1 

Si  Si  0  0 

T  -Si 

0  0  si  si 

~  ~1 

—  - 

-  _ 

(76) 


Considerable  symmetry  exists  within  the  u  matrix  as  a  result  of  the 


symmetry  of 

the 

unit 

cell 

submatrices 

and 

is  found 

"ill 

~2 

“3 

ii4 

u  = 

~2 

~1 

^4 

~3 

il4 

ill 

H2 

il4 

~3 

~2 

ill 

In  addition 

> 

iii 

=  T 

**  ‘ 

V.  T 

~i  — 

and  u  , 

~2 

and  v. 
~4 

are 

(77) 


(78) 

1,  3  the  transposes  of  the 


matrices  u. ,  v.  satisfy 
^1 


TT  =  T  u,  T 
— 1  ~  ~1  ~ 


v .  =  T  v.  T 


(79) 


These  symmetry  properties  can  be  used  to  check  intermediate  numerical  results. 
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Similarly,  the  waves  at  the  end  of  pl,te  4  c„  b  , 

°f  plate  2  via  6  re  dted  to  those  at  the  end 


-Jk-d2 


-Jk-d. 


V‘l°S  *1  ,ad  *2-  *•  defined  l„  (J„  Eo 
substituting  (au  and  <82)  £mo  “*  e'P”*"“aI  '-ora, 
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te'-A) 


in  the  form 


and 


I 

-2 

1 

=  0  . 

*1*1 

1 

-  CM 

CM 

_ i 

i  ji  are 

12x12 

matrices 

that  n 

^*3 

2  “4 

•Si 

~2 

~  ^4 

T  v3 

~2 

~1 

~2 

~2 

Zl 

V  T 

~3 i 

*1  1 

0 

0 

0 

0 

0 

0 

0 

0 

-  1 
*r 

0 

0 

0 

0 

u 

(83) 


(84) 


(85) 


where  1  is  the  3x3  unit  matrix. 

In  order  for  ©3)  to  have  non-trivial  solutions,  it  is  necessary  that 


det(u'  -  _|)  ■  0  ,  (8^) 

which  is  the  dispersion  equation  for  the  Bloc!,  waves.  Because  the  ele¬ 
ments  of  u'  depend  on  frequency  and  contains  the  wavenumbers  k^  and  k^, 

(86)  is  an  equation  of  the  form  f(u),  k  ,  k  )  ■  0,  and  can  be  solved  for  any 

x  y 

one  variable  in  terms  of  the  other  two.  In  practice,  the  equation  is  solved 
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numerically  using  a  root  search  procedure.  Numerical  aspects  of  the  solu¬ 
tions  are  discussed  below. 

B.  Numerical  Evaluation  of  the  Propagation  Characteristics 
Computations  were  made  for  a  honeycomb  material  having  Poisson's  ratio 
v  =  0.32,  which  is  the  value  for  phenolic  plastic  cores,  and  a  ratio  w/t  =  20. 
With  these  ratios  fixed,  the  normalized  wavenumbers  k  w/2  and  k  w/2  will 

x  y 

depend  on  the  shear  wave  velocity  Cg  of  the  honeycomb  material  only  through 

the  normalized  frequency  variable  kgw/ 2  =  Cgw/2uu.  At  very  low  frequencies, 

two  propagating  waves  are  found.  One  of  these  waves  has  velocity  of  propa¬ 
gation  0.862  C  that  is  independent  of  the  direction  of  propagation.  The 
s 

second  wave  has  velocity  of  propagation  that  varies  from  0.059  Cg  to  0.062  Cg 
as  a  function  of  the  direction  of  propagation.  The  5%  variation  in  the  calcu¬ 
lated  velocity  is  not  thought  to  be  meaningful,  so  that  at  low  frequencies 
the  propagation  is  isotropic  in  x-y  plane.  The  two  waves  are  referred  to  as 
the  fast  and  slow  Bloch  wave  solutions,  respectively.  For  comparison,  recall 
that  at  low  frequencies  the  Bloch  waves  with  particle  motion  along  z 

propagate  with  velocity  0.707  Cg.  For  convenience,  we  refer  to  the  Bloch  wave 

with  particle  motion  along  z  as  the  Q-SV  (quasi-shear  vertical)  wave.  The 
properties  of  the  fast  and  slow  Bloch  waves  are  discussed  separately  below. 

1.  Fast-Wave  Solution 

The  wavevector  or  slowness  curves  corresponding  to  the  fast-wave  solution 

are  depicted  in  Figure  13.  As  discussed  in  connection  with  the  Q-SV  waves,  the 

portion  of  k  -k  space  shown  in  Figure  13,  which  is  the  first  quadrant  of  the 

x  y 

first  Brillouin  zone,  contains  all  of  the  information  concerning  the  solutions 
of  (86).  The  six- fold  symmetry  of  the  dispersion  curves  about  the  z-axis  is 


Figure  13.  Slowness  curves  in  the  first  quadrant  of  the  first  Brillouin 
zone  for  the  fast  Bloch  wave  for  various  values  of  normalized 
frequency  kw/2  within  the  first  pass-band  (w/t  -  20,  v  -  0.3 


evident  from  the  plots  for  the  first  quadrant,  and  is  consistent  with  the 
hexagonal  nature  of  the  core.  The  normalized  frequencies  in  Figure  13  lie 
in  the  first  pass-band  of  the  fast-wave. 

The  dispersion  curves  for  propagation  along  x  and  along  y,  as  well  as 
their  six-fold  related  directions,  are  shown  in  Figures  14  and  15.  Above 
the  first  pass-band,  the  dispersion  curves  for  propagation  along  both  x  and 
y  consist  of  a  complex  set  of  stop  bands  and  pass-bands  that  are  not  periodic 
with  frequency,  as  in  the  case  of  SH  waves.  The  non-periodicity  results 
from  the  fact  chat  the  Bloch  wave  is  composed  of  F(0)  and  F( 1)  place  modes, 
in  addition  to  the  L  mode.  Furthermore,  the  F (0)  plate  mode  is  dispersive, 
the  coupling  at  the  V  joints  is  frequency  dependent,  and  the  interaction 
between  Y  joints  due  to  the  F(L)  mode  decreases  exponentially  with  the  square 
root  of  frequency.  The  complexity  associated  with  higher  stop  and  pass-bands 
is  discussed  in  connection  with  the  slow-wave  solution. 

As  the  frequency  increases  from  zero  to  the  first  stop  band,  the  slow¬ 
ness  curves  of  Figure  13  behave  in  a  fashion  similar  to  those  of  the  SH 
wave.  Thus,  the  curve  changes  from  being  circular  about  the  origin  to  a 
figure  that  is  slightly  elongated  along  kx>  and  six-fold  related  axes.  At 
the  normalized  frequency  ksw /2  =  0.2735,  the  slowness  curve  touches  the  bound¬ 
aries  of  the  first  Brillouin  zone  along  kx  and  related  directions.  In  the 
range  0.2735  <  kgw/2  <  0.275,  the  slowness  curve  becomes  a  series  of  closed 
shapes  about  the  apexes  of  the  Brillouin  zone. 
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While  the  fast  Bloch  wave  has  low  frequency  velocity  greater  than  the 
Q-SV  Bloch  wave,  it  has  a  lower  cutoff  frequency  for  propagation  along  x  and 
a  cutoff  band  for  propagation  along  y.  In  Figure  16  we  plotted  the  normal¬ 
ized  wavenumbers  of  the  L,  SH  and  F(0)  plate  modes,  as  well  as  the  longitu¬ 
dinal  wavenumber  of  the  honeycomb  material  versus  normalized  frequencies 
k^w/2  for  w/t  =  20.  Were  the  fast  Bloch  wave  to  be  composed  solely  of  the  L 
mode  in  the  plates,  at  low  frequencies,  one  would  expect  the  fast  Bloch  wave 
velocity  to  be  about  1.7  times  that  of  the  Q-SV  Bloch  wave;  the  value  1.7  being 
the  ratios  of  the  L  mode  to  SH  mode  velocities,  lha t  the  fast  Bloch  wave  is 
only  1.2  times  that  of  the  Q-SV  Bloch  waves  is  therefore  due  to  the  presence  of 
flexural  waves.  The  coupling  between  F(0)  and  L  modes  increases  with  frequency, 
so  that  the  phase  shift  due  to  the  very  slow  F(O)  mode  becomes  more  significant 
in  the  Bloch  wave  at  higher  frequencies.  This  is  believed  to  explain  the  rapid 
increase  in  dispersion  above  k^w/2  =  0.20,  and  a  value  of  cutoff  frequency  below 
that  of  the  Q-SV  Bloch  wave. 

In  tracing  the  higher  stop  and  pass-band  behavior  for  propagation  along  x 

it  was  found  chat  che  first  branch  of  dispersion  equation  became  complex,  as  in 

Figure  14,  and  did  not  connect  to  the  higher  pass-band  solution.  The  second 

pass-band  results  from  a  second  branch  of  the  dispersion  equation,  which  is 

imaginary  for  frequencies  below  k^w/2  =  0.341.  For  propagation  along  y,  the 

same  branch  of  the  dispersion  equation  connected  from  the  first  to  the  second 

pass-band.  The  double  solution  for  k  in  the  vicinity  of  k  w/2  *  0.6  indicates 

y  s 

that  the  slowness  curves  for  these  frequencies  consist  of  several  disconnected 
curves,  and  are  probably  quite  complex.  We  have  not  however  computed  the  shapes. 
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Tile  particle  motion  associated  with  the  Cast  Bloch  wave  is  depicted  in 
Figures  17  and  18  for  wave  propagation  along  x  and  y,  respectively,  and  a 
normalized  frequency  k  w/2  =  0.2.  These  curves  show,  in  exaggerated  scale, 
the  actual  displacement  of  points  along  the  center  line  of  the  plates  at 
some  instant  of  time.  While  the  motion  is  different  for  the  two  directions 
of  propagation,  in  each  case  it  is  symmetric  about  a  line  that  is  parallel 
to  the  direction  of  propagation  and  that  divides  the  unit  cell  in  half.  In 
addition,  the  center  of  gravity  of  the  unit  cell  appears  to  move  back  and  forth 
in  the  direction  of  wave  propagation.  Because  of  this  motion,  the  fast  Bloch 
wave  resembles  a  longitudinal  wave  in  an  elastic  continuum,  and  will  be  re¬ 
ferred  to  as  the  Q-P  (quasi- longitudina 1)  wave. 

2 .  Slow-Wave  Solution 

The  slowness  curves  corresponding  to  the  slow  Bloch  wave  are  depicted 
in  Figure  19  for  frequencies  below  the  first  stop  band  along  y.  The  cor¬ 
responding  dispersion  curves  for  propagation  along  x  and  along  y  are  shown 
in  Figures  20  and  21,  respectively.  It  is  seen  that  the  behavior  of  the 
slowness  curve  is  similar  to  that  in  Figures  8  and  13  for  the  Q-SV  and  Q-P 
Bloch  waves.  However,  the  anisotropy  effects  are  much  more  pronounced  and 
the  normalized  frequency  kgw/2  for  corresponding  curves  is  an  order  of  magni¬ 
tude  lower. 

The  difference  in  frequency  can  be  understood  by  referring  to  the  dif¬ 
ferences  in  the  wavenumbers  of  the  various  plate  modes,  as  shown  in  Figure  16. 
Over  the  ordinate  range  0  <  k ^  w/2  <  1.0,  the  flexural  wavenumber  is  more  than 
20  times  larger  than  the  shear  wavenumber.  If  the  slow  Bloch  wave  were  composed 
solely  of  the  F(0)  mode  in  the  plates,  one  would  expect  its  slowness  curve  at 
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propagation  along  y  of  the  slow  Bloch 


frequency  to  have  a  shape  similar  to  tha-  of  the  Q-SV  Bloch  wave  at  frequency 
greater  than  20  f^.  However,  the  frequency  ratio  is  somewhat  less  than  20 
due  to  the  presence  of  the  L  and  F(l)  modes. 

As  the  frequency  increases  from  zero  to  the  first  stop  band  for  propa¬ 
gation  along  x,  the  shape  of  the  slowness  curves  of  Figure  19  changes  from 
circular  about  the  origin  to  a  figure  that  is  elongated  along  k  ,  and  six-fold 
related  axes.  At  the  normalized  frequency  kgw/2  =  0.0316,  the  slowness  curve 
touches  the  boundary  of  the  Brillouin  zone.  In  the  range  0.0316  <  k^w/2  <  0.0671, 
the  slowness  curve  breaks  into  a  series  of  closed  curves  about  the  apexes  of 
the  first  Brillouin  zone.  These  curves  shrink  to  a  point  at  k^w/2  -  0.0671. 

For  frequencies  in  the  range  0.0671  <  kgw/2  <  0.112,  the  Bloch  wave  is 
evanescent  for  all  directions  of  propagation  in  the  x-y  plane. 

The  particle  motion  associated  with  the  slow  Bloch  wave  in  the  first 
pass-band  is  depicted  in  Figures  22  and  23  for  propagation  along  x  and  y 
respectively  at  normalized  frequency  kgw/2  =  0.03.  In  these  figures  we 
have  plotted,  at  an  exaggerated  scale,  the  actual  displacement  of  the 
centers  of  plates  comprising  a  unit  cell  at  a  particular  instant  in  time. 

From  Figure  22,  it  is  seen  for  propagation  along  x  the  motion  is  almost  trans¬ 
verse  to  the  direction  of  propagation.  From  Figure  23  it  is  seen  that  the 
motion  for  propagation  along  y  is  anti-symmetric  about  a  plane  through  the 
center  of  the  cell  and  parallel  to  y,  and  has  a  strong  shear  component.  On 
the  basis  of  cell  motion,  we  can  identify  the  slow  Bloch  wave  at  low  fre¬ 
quency  as  a  shear  wave  in  an  elastic  continuum,  and  will  refer  to  it  as  the 

Q-SH  (quasi-shear  horizontal)  wave. 

For  frequencies  above  the  first  pass-band,  the  dispersion  curves  of 
Figures  20  and  21  exhibit  a  complex  series  of  stop  bands  and  pass-bands,  which 
are  not  periodic  in  frequency.  For  the  case  of  propagation  along  x,  a  single 
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Figure  23.  Displacement  of  the  plates  of  a  unit  cell  at  an  instant  of 
time  during  the  propagation  of  a  slow  Bloch  wave  in  the  y 
direction  (w/t  =  20,  v  =  0.32). 


solution  to  the  dispersion  equation,  which  can  be  traced  continuously  as  a 
function  of  frequency,  gives  the  wavenumber  kx  of  propagating  waves.  How¬ 
ever,  for  propagation  along  y  it  is  necessary  to  trace  two  solutions  to  the 
dispersion  relation.  The  first  solution  gives  the  propagating  wavenumber  ky 
in  the  first  pass-band  0  <  kgw/2  <  0.0671,  while  the  second  solution  is  cutoff 

there.  At  k  w/2  =  0.112,  the  second  solution  switches  from  cutoff  to  propa- 
s 

gating,  i.e.  k  becomes  real,  and  remains  so  until  kgw/2  =  0.185.  However, 

at  k  w/2  =  0.178,  the  first  solution  again  becomes  real,  but  corresponds  to 
s 

a  solution  in  a  Brillouin  zone  above  the  first  zone  in  the  k  -k  plane.  At 

x  y 

k  w/2  =  0.185,  the  two  solutions  are  real  and  equal.  For  k  w/2  above  0.185, 
s  s 

the  two  solutions  diverge,  one  corresponding  to  the  higher  Brillouin  zone 

and  the  other  to  the  first  Brillouin  zone. 

The  foregoing  dispersion  characteristics  can  best  be  understood  by 

examining  the  slowness  curves  as  drawn  for  a  progression  of  frequencies  about 

k  w/2  *  0.185.  As  k  w/2  increases  from  k  w/2  =  0.112,  the  slowness  curve 
s  s  s 

starts  as  a  circle,  as  shown  in  Figure  24,  whose  radius  increases  with  fre¬ 
quency.  The  shape  of  the  circle  distorts  with  increasing  frequency  until  the 

slowness  curve  touches  the  boundary  of  the  first  Brillouin  zone  at  k  w/2  = 

s 

0.178,  as  shown  in  Figure  24.  Because  the  slowness  curve  touches  the  boundary 

at  intermediate  angles,  rather  than  along  the  k^  and  related  six-fold  axes, 

for  frequencies  above  0.178  it  takes  the  form  of  closed  curves  about  the 

vertices  of  the  Brillouin  zone  and  about  the  points  k  *  rr/3,  k  =0  and  six- 

x  y 

fold  related  points.  In  the  range  0.178  <  kgw/2  <  0.185,  the  closed  curves 

about  the  vertices  shrink  to  a  point,  and  then  expand  above  0.185.  The  closed 

curves  about  k  •  tt/3 ,  k  *  0,  and  related  points,  decrease  in  size  in  the 
x  y 
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range  0.178  <  k  w/2  <  0.188.  For  k  w/2  between  0.188  and  0.207  no  real  solu- 
s  s 

tion  for  k^  exist,  as  is  seen  from  Figure  20. 

The  closed  curves  about  the  vertices  of  the  Brillouin  zone  expand  in  size 

for  k  w/2  >  0.185,  as  seen  in  Figure  25,  until  they  touch  the  k  and  related 
s  x 

axes  at  kgw/2  =  0.207.  Above  0.207,  the  slowness  curve  is  again  a  closed 
curve  about  the  origin,  whose  radius  decreases  with  frequency  to  zero  at 
kgw/2  =  0.0342.  The  foregoing  discussion  illustrates  the  complex  dependence 
of  the  Bloch  wave  properties  on  frequency  in  the  higher  pass-bands. 
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V.  BLOCH  WAVE  PROPAGATION  OBLIQUE  TO  THE  CELL  AXES 


For  propagation  oblique  to  the  cell  axes,  the  scattering  matrix 
for  the  Y  joint  is  1  2x1  2  and  the  matrices  au  and  lx  are  1x4,  as  de¬ 
fined  in  (38).  vVith  this  recognition,  the  matrix  algebra  needed  to 
find  the  propagation  characteristics  of  the  Bloch  wave  is  the  same  as 
that  developed  in  Section  IV,  provided  the  matrices  defined  there  are 
appropriately  generalized.  Thus,  the  matrices  S^  in  (58)  that  are 
used  to  partition  S  are  4x4,  as  are  the  matrices  S.'  in  (61).  The  e 
matrix  defined  in  (60)  must  be  generalized  as 


where  Kg,  Kq  and  K  j  are  the  modal  wave  numbers  along  x,  as 

defined  in  (7),  (9)t  (11)  and  (15). 

Because  the  SH  mode  is  symmetric,  it  has  the  same  junction  pro¬ 
perty  (66)  as  the  L  mode  so  that  the  T  matrix  defined  in  (7  2)  becomes 

10  0  0 

0  10  0 

(88) 

0  0-10 
0  0  0  -1 
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Finally,  in  the  ^  matrix  of  (85)  one  must  use  the  4x4  unit  matrix  for 
1^  and  recognize  that 


U,  2=  exp 


(89) 


The  wave  number  kz  is  contained  in  the  scattering  matrix  S  and 

hence  in  u\  so  that  the  dispersion  equation  (86)  is  of  the  form 

f(o>.  k^,  k^)  =0.  In  the  rest  of  this  section  we  report  on  numerical 

solutions  of  the  dispersion  relation  for  k  vs.  k  with  k  =  0  and  fixed 

x  z  y 

values  of  w,  and  for  k  vs.  k  with  k  =0  and  fixed  values  for  u>. 

y  z  x 

These  results,  together  with  those  of  the  previous  section,  indicate 

the  shape  of  the  slowness  surface  by  means  of  cuts  in  the  (k  ,k  ) 

’  x  y 

plane,  in  the  (k  ,  k  )  plane  and  six -fold  related  planes,  and  in  the 
X  z 

(k  ,  k  )  and  six -fold  related  planes.  Calculations  were  carried  out  for 
y  z 

a  honeycomb  material  having  Poisson's  ratio  v  =  0.  3  2  and  for  cell- 
wall  aspect  ratio  w/t=  20. 


A.  Low  Frequency  Propagation:  Continuum  Model 

For  normalized  frequencies  below  about  kgw,/ 2  =  0.  01,  all  of  the 

Bloch  waves  propagating  in  the  (kx>  k^)  plane  show  linear  variation 

of  k  and  k  with  w,  and  hence  are  non-disper sive.  We  have  there- 
x  y 

fore  computed  the  slowness  curves  for  kgw/ 2  =  0.  01,  and  plotted 
kz/kg  versus  or  k^,/kg.  The  resulting  curves  are  independent 

of  frequency  for  kgw/2  <  0.  01. 

In  Figure  26  we  have  plotted  the  slowness  curves  for  the  two 
fastest  Bloch  waves.  The  same  numerical  values  were  obtained  for 
propagation  in  the  (x,  z)  plane  as  for  propagation  in  the  (y,  z)  plane. 
The  inner  curve  is  a  quasi-longitudinal  wave  (Q-P),  as  noted  from 
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the  particle  displacement  for  propagation  along  x  and  y.  Also,  for 

propagation  along  z,  the  value  k^  kg  =  0,  59  is  the  same  as  for  the 

wavenumber  kT  of  the  symmetric  Lamb  mode  in  the  cell  walls,  which 

is  essentially  an  extensional  or  longitudinal  wave  at  low  frequencies. 

Based  on  the  polarization  properties  for  propagation  in  the  (x,  y)  plane, 

the  outer  curve  in  Figure  26  represents  a  vertically  polarized  shear 

wave  { Q-SV).  The  fact  that  the  intersection  of  the  curve  with  the 

k  /k  or  k  /k  axis  and  with  the  k  /k  axis  have  the  same  values 
x  s  y  s  z  s 

(41)  indicates  that  the  honeycomb  acts  as  an  elastic  continuum^ 

Because  the  slow  Bloch  wave  travels  at  one  tenth  the  velocity  of 
the  fast  Bloch  waves  in  the  (x,  y)  plane,  its  slowness  curve  cannot  be 
drawn  with  equal  scales  for  k^  and  k^,  as  in  Figure  26.  For  this 
reason  we  have  drawn  the  dispersion  curves  with  different  scales  in 
Figure  27.  The  slow  Bloch  wave  is  labeled  Q-SH  because  its  parti¬ 
cle  motion  is  similar  to  horizontally  polarized  shear  wave  for  propa¬ 
gation  along  x  or  y.  The  discrepancy  between  the  slowness  curves 
for  propagation  in  the  (x,  z)  planes  and  (y,  z)  plane  amounts  to  2.  5% 
to'-,  is  thought  to  result  from  the  approximations  used  to  compute 

•  it  is  seen  that  both  shear  waves  cross  the  k  axis  at  the 

z 

.  ■  nt,  in  accordance  with  the  properties  of  a  hexagonal  'ntin- 

■  •  •'  examine  further  the  correspondence  between  the  slow- 
Kigures  26  and  27  and  an  elastic  continuum  of  hexa- 
v <•  have  used  the  propagation  characteristics  along 
-  •  define  the  stiffness  constants  of  the  continuum. 

•  »  'he  slowness  curves  for  continuum  were  com- 
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puted  and  compared  with  those  of  Figures  26  and  27. 

Static  measurements  of  the  mass  density  and  stiffness  con- 

(12) 

stant  a  WRXI  kraftboard  honeycomb'  ,  as  listed  in  Table 

III,  were  used  in  conjunction  with  the  properties  listed  in  Reference 

(11)  to  find  the  remaining  stiffness  constants,  which  are  listed  in 

Table  III.  The  mass  density  p  of  the  kraftboard  listed  in  Table  III 

was  obtained  by  dividing  by  the  fraction  of  honeycomb's  volume 

occupied  by  the  phenolic  material.  Assuming  the  kraftboard  act  as  an 

isotropic  medium  with  v  =  0.  32,  and  recognizing  that  k  for  the  Q-P 

z 

wave  propagating  along  z  is  equal  to  k^,  (j  can  be  found  from  the  re¬ 
lation 


Knowing  ^  and  v  ,  one  can  find  X.  from  the  expression 

-  V 

V  "  2(\+fx) 

for  Poisson's  ratio.  These  values  are  also  listed  in  Table  III. 

The  average  of  the  values  of  kx/kg  was  used  in  obtaining  Cgg 
in  Table  III.  With  the  values  listed  in  Table  III,  we  have  computed 
the  slowness  curves  for  a  hexagonal  continuum  and  compared  them 
with  Figures  26  and  27.  The  Q-P  and  Q-SV  curves  were  indistinguish¬ 
able,  while  the  Q-SH  curve  fell  between  the  k  vs.  k  and  k  vs.  k 

z  x  z  y 

curves  of  Figure  27,  We  therefore  conclude  that  the  continuum  model 
is  an  accurate  representation  for  the  honeycomb  at  low  frequencies. 

It  should  be  noted  that  for  a  WR11  kraftboard  honeycomb,  the  nor¬ 
malized  frequency  kgw/2=  0.01  corresponds  to  an  actual  frequency 
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TABLE  III 


material  constants  for  a  wrii  kraftboard  honeycomb 


Kraftboard 

Honeycomb 

p  =  0.  899  gm/ cm5 

=  0.061  gm/cm5 

v  =  0.  3  2 

=  3.  79x10^  dynes/cm2 

H  =  3.  38x1 0i0  dyne/ cm 2 

Ct  =  9.  75xl08 

\  =  1. 9x1 010 

Cl2=  9.66X108 

C  3  =  5.  94x1 0S 

C44  =  6-  60x1 08 

C„„  =  4.  60x1 06 

66 
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f  —  '30  Hz.  For  higher  frequencies,  the  Q-SH  wave  becomes  disper¬ 
sive,  while  the  Q-P  and  Q-SV  waves  do  not  exhibit  dispersive  effects 
until  a  frequency  of  nearly  10  KHz. 


B.  Dispersive  Effects  for  the  Slow  Bloch  Wave 
To  illustrate  the  dispersive  effects  on  the  slow  Bloch  wave  we 
have  computed  slowness  curves  for  kgw/ 2  =  0.03,  0.06,  0.13  and  0.26. 
The  results  for  kgw/2=  0.03  are  plotted  in  Figure  28.  The  extent 
of  the  curves  along  and  is  markedly  different.  In  the  insert 
of  Figure  28  we  have  sketched  a  perspective  view  of  the  slowness 
surface  in  the  first  Brillouin  zone.  The  surface  is  seen  to  have  the 


shape  of  a  six-pointed  star.  Slowness  surfaces  in  the  neighboring 
Brillouin  zones  have  the  same  shape,  and  oriented  such  that  the 
points  of  the  stars  face  those  of  the  star  in  the  first  zone. 

As  the  frequency  increases  from  0.0  3,  the  tips  of  the  stars  in 
neighboring  Brillouin  zones  grow  towards  each  other  until  they  touch 
at  0.316  and  form  tunnels  for  frequencies  above  0.316.  The  tunnel  shape 
is  indicated  in  Figure  29,  which  is  drawn  for  kBw/2  =  0.06.  As  the 
frequency  increases  further,  the  tunnels  increase  in  width  until  they 
merge  at  0.071  when  the  wave  becomes  cut-off  along  y,  as  indicated  in 
Figure  21.  For  frequencies  between  0.071  and  0.112,  where  the  second 


pass  band  starts,  the  slowness  surface  in  (k^,  k^,  k^)  space  undulates 

but  does  not  touch  down  to  the  (k  ,  k  )  plane. 

x  y 

The  slowness  surfaces  of  the  slow  Bloch  wave  in  the  second  pass 


band  consist  of  two  surfaces.  The  outer  surface  does  not  touch  the 


(k  ,  k  ),  while  the  inner  surface  in  the  first  Brillouin  zone  is  closed 
x  y 

about  the  origin.  Identical  closed  surfaces  exist  in  each  of  the 
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Brillouin  zones.  The  intersections  of  these  surfaces  with  the  (k  .  k  ) 

x  z 


and  (k  ,  k  )  planes  are  shown  in  Figures  30  and  3l  for  k  w/  2  =  0.  1  3. 
y  z  s 

The  outer  and  inner  surfaces  touch  in  the  (k  ,  k  )  plane  and  six-fold 

y  Z 

related  planes.  The  touching  points  are  like  the  points  of  conical 
refraction  in  the  optical  slowness  surfaces  of  a  biaxial  crystal.  The 
surfaces  for  kz  >  0  may  be  thought  of  as  a  head  (inner  surface)  with  a 
hat  (outer  surface)  that  rests  on  the  head  at  six  points. 

As  an  example  of  the  slowness  curves  in  the  third  pass  band  of 
the  SH-Bloch  wave,  we  have  drawn  Figures  3  2  and  3  3  for  kgw/2  =  0.  26. 
Slowness  curves  for  propagation  in  the  (x,  z)  and  six -fold  related  planes 
are  indicated  in  Figure  3  2,  while  Figure  33  has  been  drawn  for  propa¬ 
gation  in  the  (y,  z)  and  six-fold  related  planes.  The  Q-SH  slowness 
surface  is  seen  to  consist  of  two  parts.  The  upper  part  represents 

an  undulating  surface  that  does  not  touch  the  (k  ,  k  )  plane.  The  lower 

x  y 

Q-SH  slowness  surface  resembles  a  tube  enclosing  the  boundary  of 
the  Brillouin  zone.  The  character  of  the  slowness  curves  for  the  fast 
Bloch  waves  is  discussed  in  the  next  section. 


C.  Dispersive  Effects  for  the  Fast  Bloch  Wave 
As  seen  from  Figures  9  and  1  0,  the  Q-SV  Bloch  wave  propagat¬ 
ing  in  the  (x,  y)  plane  does  not  have  significant  dispersion  for 
kgw/  2  <  0.  4.  The  Q-P  Bloch  wave  propagating  in  the  (x,  y)  plane  how¬ 
ever  shows  significant  dispersion  even  at  kgw/2  =  0.  2.  For  oblique 
propagation  the  dispersive  effects  should  be  evident  in  both  waves, 
even  at  frequencies  as  low  as  kgw/2  =  0.  2,  since  both  waves,  are 
composed  of  the  same  plate  modes  in  the  honeycomb  walls. 
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k  /k 
z  s 


Figure  30.  Slowness  curves  in  the  (kx,kz)  plane  for 

frequency  kgw/2  ■  0.13  in  the  second  p*BS- 
band  of  the  Q-SH  Bloch  wave  showing  the 
presence  of  a  fourth  surface  (w/ t  ■  20, 
v  -  0.32). 
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1.22  1.42  3.61 

Figure  31.  Slowness  curves  in  the  (k  ,k  )  plane  for  frequency 
ksw/2  ■  0.13  in  the  second  piss-band  of  the  Q-SH 
Bloch  wave  showing  the  presence  of  a  fourch  surface 
(w/ t  =  20,  v  ■  0.32).  Ihis  figure  is  for  the  same 
frequency  as  Figure  30. 
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Figure  32.  Slowness  curves  in  the  (kx,kz) 
k  w/2  =  0.26  in  the  third  pass- 
Bioch  wave  (w/ t  *  20,  v  =  0.32) 


4.03 

plane  for  frequency 
band  of  the  Q-SH 


The  nature  of  the  dispersive  effect  can  be  understood  by  exam¬ 
ining  propagation  for  frequencies  in  the  vicinity  of  kgw  2  =  0.  21  6, 
where  the  wavenumbers  of  the  Q-P  and  Q-SV  Bloch  waves  propagating 
in  the  (x,  y)  plane  are  equal.  The  slowness  curves  of  the  Q-P  and 

Q-SV  Bloch  waves  are  drawn  in  Figure  34  for  k  w  ■  2  =  0.  21.  Because 

s 

anisotropy  effects  are  not  significant  at  this  frequency,  the  same 

curve  holds  for  propagation  in  the  (x,  z)  plane  and  for  propagation  in 

the  (y,  z)  planes.  Comparing  Figure  34  with  Figure  26,  which  is 

drawn  for  k  *w  /  2  =  0.01.  it  is  seen  that  the  Q-P  wave  shows  disper- 
s 

sive  effects  for  propagation  in  the  (x,  y)  plane.  However,  for  propa¬ 
gation  along  z,  the  Q-P  wave  shows  no  dispersion  since  its  wave- 
number  is  that  of  the  symmetric  Lamb  wave  in  the  cell  wall.  For 
this  direction  of  propagation  it  is  the  Q-SV  wave  that  exhibits  disper¬ 
sive  effects. 

At  the  frequency  kgw/2=  0.  23,  the  slowness  curves  shown  in 
Figure  35  reflect  the  fact  the  wavenumber  of  the  Q-P  wave  propagat¬ 
ing  in  the  (x,  y)  plane  is  greater  than  that  of  the  Q-SV  wave.  However, 
for  propagation  along  z,  the  relative  size  of  the  wavenumbers  is  re¬ 
versed.  Thus,  along  the  inner  branch  of  the  slowness  curve  in  Figure 
35  the  Bloch  wave  polarization  switches  from  that  of  a  shear  wave 
propagating  in  the  (x,  y)  plane  with  particle  motion  along  z  to  that  of  a 
longitudinal  wave  propagating  along  z.  Alternatively,  the  outer  branch 
has  the  polarization  of  a  longitudinal  wave  for  propagation  in  the  (x,  y) 
plane,  since  then  the  SH  plate  modes  are  not  part  of  the  Bloch  wave. 
However,  for  propagation  along  z,  the  outer  branch  has  shear 
polarization. 


Thus,  in  the  first  pass  band  of  the  fast  Bloch  waves,  a  significant 
aspect  of  the  dispersion  is  to  cause  the  polarizat  ion  type  of  the  wave 
to  change  with  direction  of  propagation.  As  a  result,  the  Bloch  waves 
cannot  be  labeled  as  a  quasi-shear  or  quasi-longitudinal  for  frequencies 
above  kgw/2=  0.  216.  This  conclusion  is  further  supported  by  the 
slowness  curves  for  the  fast  Bloch  waves  at  the  frequency  kgw/2=  0.  26, 
which  are  shown  in  Figures  3  2  and  33.  From  Figures  30  and  31  it  is 
seen  that  the  dispersion  of  the  fast  Bloch  waves  first  becomes  significant 
at  kgw/2  =  0.  13  for  shear  wave  progagation  along  z. 

Because  the  fast  Bloch  waves  do  not  show  significant  dispersion 
for  normalized  frequencies  below  kgw/ 2  =  0.  1  3,  their  propagation 
can  be  modeled  by  that  of  Q-P  and  Q-SV  plane  waves  in  a  hexagonal 
continuum  for  kgw/2  <  0.  1  3.  However,  the  continuum  model  can  be 
used  for  Q-SH  wave  propagation  only  in  the  range  kgw/  2  <  0.  013. 

Thus,  in  dealing  with  bounded  honeycombs,  the  continuum  model  can 
be  used  for  k  w/2  <  0.  013  unless  Q-SH  waves  are  not  excited,  in 
which  case  the  continuum  model  can  be  used  up  to  kgw/  2  =  0.  13.  This 
distinction  is  used  in  the  next  section. 
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VI.  WAVES  GUIDED  BY  A  FACE  PLATE  BONDED  TO  A 
HONEYCOMB 


A  complete  treatment  of  wave  propagation  in  honeycomb  panels 
requires  taking  into  account  both  of  the  face  plates  and  the  finite 
honeycomb  thickness.  The  properties  of  these  waves  result  from  a 
combination  of  effects  due  to  coupling  between  the  face  plates  and 
the  honeycomb,  and  to  the  thickness  of  the  honeycomb.  However, 
for  the  NDE  task  of  finding  areas  of  debonding,  one  is  only  interested 
in  coupling  between  the  face  plate  and  honeycomb.  This  coupling  can 
be  studied  in  isolation  by  examining  the  loading  that  a  semi-infinite 
honeycomb  places  on  the  modes  of  the  face  plate. 

In  this  section  we  consider  the  properties  of  the  lowest  modes 
guided  by  an  aluminum  face  plate  that  is  bonded  to  a  semi- infinite  honey 
comb,  as  shown  in  Figure  36.  The  honeycomb  is  assumed  to  act  as 
a  hexagonal  continuum  having  the  parameters  listed  in  Table  III  for  a 
Shelter  Core  honeycomb.  The  face  plate  is  assumed  to  have  thickness 
h  =  1  / 1  6  in.  and  to  be  bonded  to  the  honeycomb  by  means  of  an  ad¬ 
hesive  layer  of  vanishing  thickness.  However,  in  practical  cases  the 
adhesive  layer  can  have  major  influence  on  the  guided  wave,  as  is 
shown  in  the  next  section  for  foam-core  panels. 

The  SH  mode  in  the  aluminum  plate  couples  to  the  slow  Q-SH 
Bloch  wave  in  the  honeycomb.  However,  both  the  symmetric  Lamb 
mode  and  the  flexural  mode  in  the  aluminum  plate  couple  to  the  fast 
Bloch  waves.  Thus  we  may  treat  the  SH  modes  separately  from  the 
Lamb  and  flexural  modes. 

The  Q-SH  Bloch  wave  is  non-dispersive  for  k8w/2  less  than  about 
0*  013,  using  the  dimensions  assumed  for  the  Shelter  Core  honeycomb 
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gives  an  upper  frequency  limit  of  about  1,  000  Hz.  Thus  our  study 
for  guided  SH  waves  in  the  aluminum  plate  is  limited  to  frequencies 
less  than  1,  000  Hz.  On  the  other  hand,  the  fast  Bloch  waves  are 
non-dispersive  for  kgw/2  <0.  13  or  f<  10  kHz,  so  that  the  continuum 
model  for  the  fast  Bloch  waves  may  be  used  up  to  this  frequency. 
Hence,  our  study  of  the  symmetric  Lamb  and  flexural  modes  of  the 
aluminum  plate  apply  up  to  1  0  kHz. 


A.  SH  Mode  of  the  Aluminum  Plate 

The  fields  associated  with  the  SH  mode  in  the  aluminum  plate  can 
be  represented  as  a  combination  of  the  SH  plane  waves  in  aluminum 
propagating  with  wave  vectors  (k  +  k1  )  where 


k' 


zs 


/<*)2-kx2 


(90) 


and  k*  is  the  shear  wavenumber  of  aluminum.  The  wavenumber  k 
s  x 


along  x  must  be  the  same  for  both  plane  waves,  and  for  the  Q-SH  Bloch 
wave  in  the  honeycomb  in  order  to  satisfy  the  boundary  conditions  at 
the  free  surface  and  at  the  plate -honeycomb  boundary.  The  field  quan¬ 
tities  entering  into  the  boundary  conditions  at  z  =  +  h/  2  are  v  and 
T  .  If  A+  and  A  are  the  amplitudes  of  these  two  waves,  then  the 
field  quantities  in  the  plate  are  given  by  (3)s 


+  "jkzsz  -  jkzsZ  ”jkxx 
V  (x.  *)  =  (Ae  +  A  e  ZS  )  e  x 


and 


T  (x, 
yz' 


Ll’k’  ,  -jk*  z  ik*  z  -jk  x 

*  z>  =  —  (-A  e  +  A  e  )  e 

G) 


(91) 


(92) 


where  is  the  shear  modulus  of  aluminum. 


If  D"  represents  the  amplitude  of  the  Q-SH  wave  in  the  honeycomb, 

(13) 

then  the  field  quantities  for  z  <  -h/2  are  ' 


jk  z  -jk  x 

vy(x,  z)  =  D“e  z  e  x 


and 


C44kz  -  jkzz  "ikxX 
T  (x,  z)  s  -  -  Dee  x 

yZ  '  ^  Cl) 


(93) 


(94) 


Here  C,.  is  the  stiffness  constant  listed  in  Table  III  and  k  is  the 
44  z 

wavenumber  along  z  of  the  Q-SH  wave.  Using  the  continuum  model  for 
the  honeycomb,  kz  can  be  found  from  the  dispersion  equation 


C66kx  +  C4A2-'>H“2=  0  (95) 

where  and  Cj.g  are  taken  from  Table  III. 

The  boundary  condition  at  the  free  surface  z  =  h/2  is  Tyz(x,  h/2)  =  0, 
while  at  the  interface  z  =  -h/2  we  have  that  v  (x,  z)  and  T  (x,  z)  are 

y  yz 

continuous.  After  some  manipulation,  these  conditions  can  be  written 
in  the  form 


> 

+ 

_ i 

iHsh  !  ,  i  2  K.h'2 

e  i  1  i-  e  e 

A" 

•iki,h  j  :  c66k,  2  )k,«h/2 

-e  J  1  J-  ,  ,  «  e 

D 

J  i  *k,s 

L 

(96) 


The  dispersion  equation  for  the  wave  is  obtained  by  setting  the  deter¬ 
minant  of  the  3x3  matrix  equal  to  zero.  Recognizing  that  kz  and  k’zs 
are  functions  of  k^,  the  dispersion  equation  takes  the  form  f(kx,«)  =  0. 
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Because  the  Q-SH  Bloch  wave  is  extremely  slow,  for  kx  in  the 

i 

vicinity  of  k  ,  this  Bloch  wave  will  be  propagating  along  z.  As  a 
s 

result,  the  SH  mode  in  the  aluminum  plate  will  be  converted  into  a 
leaky  surface  wave  that  radiates  its  energy  into  the  honeycomb.  The 
continuous  radiation  of  energy  causes  an  exponential  decrease  in  the 
amplitude  of  the  surface  wave  with  distance  in  the  propagation  direc¬ 
tion.  The  exponential  decrease  exhibits  itself  as  an  imaginary  part 
in  the  wavenumber  k^  of  the  surface  wave.  Thus  the  root  of  the  dis¬ 
persion  equation  f(kx<o>)  =  0  is  the  complex  number  kx  =  (3  -  jo ,  where 
3  is  the  phase  constant  and  a  the  attenuation  constant  of  the  leaky 
surface  wave. 

Numerical  solutions  of  the  dispersion  equation  for  3  and  a  are 
plotted  in  Figure  37.  The  leaky  wave  is  somewhat  slower  than  the  SH 
wave  in  an  isolated  aluminum  plate.  Because  a  and  3  are  on  the  same 
order,  the  wave  is  heavily  attenuated  even  in  propagating  one  wave¬ 
length.  Thus,  the  coupling  to  the  honeycomb  strongly  affects  the  wave 
in  the  aluminum  plate.  In  principle,  the  significant  perturbation  re¬ 
sulting  when  the  plate  is  bonded  to  the  honeycomb  makes  possible  detec¬ 
tion  of  debonded  regions.  However,  the  wavelength  at  1  KHz  is  2.  4 
meters,  which  indicates  that  only  large  areas  could  be  inspected.  To 
achieve  high  resolution,  it  is  necessary  to  use  much  higher  frequency 
waves,  for  which  the  continuum  model  of  the  honeycomb  is  not  valid. 

B,  Lamb  Modes  of  the  Aluminum  Plate 

When  treating  the  symmetric  and  anti -symmetric  (flexural)  Lamb 
modes,  it  is  necessary  to  include  SV  and  P  plane  waves  in  the  plate 
that  have  components  of  propagation  in  the  +z  direction  as  well  as  in 
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,  a  (m  ) 


Attenuation  constant  cr  and  propagation  constant  0  of 
the  SH  leaky  wave  guided  by  an  aluminum  face  plate 
bonded  to  a  semi-infinite  honeycomb. 


the  -z  direction.  The  fields  in  the  honeycomb  are  composed  of  the 

Q-P  and  Q-SV  Bloch  waves  propagating  in  the  -z  direction.  The  am- 

+  + 

plitude  of  these  six  waves  are  designated  as  A  ,  B  ,  C  and  D  in 
Figure  36.  The  boundary  conditions  on  v^,  v^,  T  and  T zz  at 
z  =  +  h  /  2  requires  that  all  of  the  plane  waves  have  the  same  wave- 
number  kx  along  x.  By  requiring  that  the  plane  waves  satisfy  the 
boundary  conditions,  one  obtains  the  dispersion  equation  for  the  guided 
waves. 


1.  Formal  Solution  for  the  Guided  Waves 


The  fields  in  the  plate  are  given  by  the  sum  of  four  plane  wave 
constituents.  The  particle  velocity  is  given  by^ 
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where  kg  and  kp  are  the  wavenumbers  of  the  shear  and  longitudinal 
waves  in  the  aluminum.  The  components  of  stress  entering  into  the 
boundary  conditions  are 


T  (x,  z)  = 
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and 
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where  X.  is  the  Lame  coefficient  of  the  aluminum. 

To  obtain  the  fields  in  the  honeycomb,  we  make  use  of  the  pro¬ 
perties  of  a  hexagonal  continuum^  K  Thus,  the  wavenumbers  k 

zs 

and  kzp  of  the  Q_s  and  Q-P  waves  can  be  found  from  the  dispersion 
equation 


(Cllk*  +  C44kz  '  pH“2)<C44k*  *  C33kz  ‘  ^h"2* 

-(Cl3  +  C44>2k,\2-°- 


The  components  of  particle  velocity  are 
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and 
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Finally,  the  components  of  stress  entering  the  boundary  conditions  can 
be  found  from 
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and 
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The  boundary  conditions  at  z  =  h/ 2  require  that  T  (x,  h/ 2)  and 

xz 

Tzz(x.  h/2)  vanish.  At  the  interface  z  =  -h/ 2,  the  velocity  components 

v  (x,  -h/2)  and  v  (x,  -h/2)  together  with  the  stress  components 
X  z 

T  (x,  -h/2)  and  T  (x,  -h/2)  must  be  continuous.  These  conditions 
xz  zz 

take  the  form  of  six  homogeneous  equations  in  the  six  unknown 
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amplitude  coefficients.  These  equations  can  be  written  in  the  form 
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where  the  6x6  matrix  L  is  given  in  Table  IV.  The  elements  in  the 
matrix  depend  on  u  and  kx.  Hence,  setting  the  determinant  of  L 
equal  to  zero  in  order  for  (109)  to  have  non-trivial  solutions,  gives 
the  dispersion  equation  for  the  guided  waves  in  the  form  f(kx,  w)  =  0. 

2.  Numerical  Evaluation  of  the  Dispersion  Characteristics 

Results  obtained  for  the  dispersion  characteristics  of  the  leaky 
Lamb  mode  are  summarized  in  Figure  38,  where  we  have  plotted  the 
phase  constant  |3  and  attenuation  constant  a  as  a  function  of  frequency. 
For  comparison  we  have  also  plotted  the  wavenumber  k' T  of  the  sym- 
metric  Lamb  mode  of  a  free  aluminum  plate.  It  is  seen  that  the  honey¬ 
comb  slows  the  wave  so  that  0  >  k'^  and  causes  it  to  radiate  or  leak 
into  the  honeycomb.  The  attenuation  becomes  constant  with  frequency 
above  2KHz  and  corresponds  to  a  l/e  attenuation  length  of  0.8  m.  At 
lOKHz,  the  wavelength  Zw/fi  =  0.56  m,  so  that  wave  is  strongly  attenu* 
ated  per  wavelength. 
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ELEMENTS  OF  THE  MATRIX  L 


1 


Attenuation  constant  a  and  propagation  constant  &  of  the  leaky  Lamb  wave  guided 
by  an  aluminum  face  plate  bonded  to  a  semi- infinite  honeycomb. 


The  flexural  mode  of  the  aluminum  plate  is  slow  compared  to  the 
Q-SH  and  Q-SV  waves  of  the  honeycomb  for  frequencies  below  10  KHz.  As 
a  result,  it  is  expected  that  this  mode  will  not  leak  into  the  honeycomb, 
but  remain  a  bound  surface  wave  in  the  plate.  Ihis  supposition  is  borne 
out  by  the  values  obtained  for  the  root  of  the  dispersion  equation,  from 
which  it  is  found  that  the  attenuation  constant  a  =  0.  In  Figure  39  we 
have  plotted  g  obtained  from  the  dispersion  equation.  For  comparison  we 
have  plotted  the  wavenumber  k^  of  the  flexural  wave  in  a  free  plate  and 
the  wavenumber  k^  of  the  Rayleigh  wave  propagating  on  a  z-cut  hexagonal 
continuum  having  the  elastic  constants  taken  for  the  honeycomb  at  low  fre¬ 
quencies.  Below  3  KHz  the  surface  wave  is  essentially  a  Rayleigh  wave  in 
the  honeycomb.  At  higher  frequencies  the  plate  acts  as  a  perturbation  on 
the  honeycomb  and  0  deviates  from  k^.  Ihe  surface  wave  is  slow  compared 
to  the  leaky  Lamb  mode  shown  in  Figure  38.  At  10  KHz,  the  wavelength  of 
the  surface  wave  is  0.057  m  or  5.7  cm,  which  is  an  order  of  magnitude 
smaller  than  that  of  the  Iamb  mode. 

For  application  to  NDE,  it  is  necessary  that  the  wavelength  be  on 
the  order  of  the  cell  diameter,  which  is  about  1  cm.  Thus,  the  flexural 
mode  would  be  useful  at  frequencies  somewhat  above  10  KHz.  However,  the 
Lamb  and  SH  modes  would  require  frequencies  above  100  KHz  to  achieve  such 
small  wavelengths.  Such  high  frequencies  lie  in  higher  pass  or  stop  bands 
of  the  honeycomb,  where  the  periodicity  of  the  honeycomb  is  significant, 
and  must  be  accounted  for  in  computing  the  leaky  wave  characteristics.  It 
therefore  seems  most  convenient  to  use  the  flexural  mode  for  NDE  applica¬ 
tions.  This  mode  also  offers  the  advantage  of  a  large  difference  between 
and  0,  so  that  debonded  regions  of  the  honeycomb  are  easily  distinguished 
from  bonded  regions. 
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VII.  FOAM  CORE  PANELS 

We  consider  here  wave  propagation  in  form  core  panels  in  the 
frequency  rang  100  KH  to  1  MH  for  application  to  NDE.  In  this 
frequency  range,  the  acoustic  attenuation  in  the  foam  is  large.  As  a 
result,  elastic  waves  excited  in  one  face  plate  will  be  highly  attenuated 
by  the  time  they  reach  the  opposite  face  plate.  Because  of  this  pro¬ 
perty,  the  wave  propagation  can  accurately  be  determined  by  consider¬ 
ing  a  single  face  plate  bonded  to  a  semi-infinite  foam  block,  as  was 
done  for  the  honeycomb  in  the  previous  section.  For  this  geometry 
we  examine  the  perturbation  due  to  the  foam  in  the  modes  guided  by  the 
face  plate. 

Two  methods  exist  for  bonding  the  foam  to  the  aluminum.  For  one 

method  the  chemical  agents  used  to  produce  the  foam  are  allowed  to 

react  directly  on  the  aluminum.  In  the  process  of  reacting,  the  foam 

becomes  bonded  to  the  aluminum.  In  the  second  method,  preformed 

foam  slabs  are  glued  to  the  aluminum.  This  second  method  introduces 

an  additional  glue  layer  whose  elastic  properties  greatly  influence  the 

(14) 

characteristics  of  the  guided  wave  .  in  what  follows,  we  consider 
separately  the  aluminum-foam  case  and  the  aluminum  -  glue  -  foam 
configur  ation. 

It  is  assumed  that  the  aluminum  plate  has  thickness  l/l6  in.  and 

the  elastic  constants  used  for  the  foam  are  the  nominal  values  listed 

(151 

for  rigid  polyurethan  foams'  '.  The  glue  was  assumed  to  be  a 
(15) 

syntactic  epoxy  .  Elastic  constants  used  for  the  various  materials 
are  listed  in  Table  V.  The  method  of  analysis  used  to  find  the  propaga¬ 
tion  constants  of  the  guided  waves  is  the  same  as  that  discussed  in  the 
previous  section.  For  the  aluminum  -  foam  case,  the  expressions  of 
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the  previous  section  can  be  used  directly  by  substituting  the  mass 
density  of  the  foam  for  and  setting  =  ^33  =  *■£  +  2  = 

CJ3  =  X^.  and  where  X^  and  are  the  Lame  constants  for  the 

foam.  For  the  aluminum  -  glue  -  foam  configuration,  a  similar 
analysis  is  applied  that  accounts  for  the  plane  waves  in  the  glue  layer. 

A,  Aluminum  -  Foam  Configuration 

In  Table  V  we  have  listed  the  shear  and  longitudinal  wave  veloci¬ 
ties  of  the  materials,  as  computed  from  the  elastic  constants  that  are 
listed  there.  It  is  seen  that  the  wave  velocities  in  foam  are  much 
smaller  than  those  in  aluminum,  so  that  the  SH  and  symmetric  Lamb 
modes  guided  by  the  plate  will  radiate  or  leak  energy  into  the  foam. 
Over  the  frequency  range  considered  here  even  the  flexual  mode  of  the 
plate  is  fast  compared  to  the  waves  in  the  foam,  so  that  this  wave  will 
also  be  of  the  leaky  type. 

The  wavenumber  k^of  the  leaky  waves  takes  the  form  kx  =  P-jcr  , 
where  3  is  the  phase  constant  and  a  the  attenuation  constant  that  results 
from  the  radiation  into  the  foam.  Because  of  the  low  density  of  the 
foam,  its  presence  does  not  significantly  perturb  the  value  of  3.  In 
Figure  40  we  have  plotted  the  phase  constant  3  for  the  lowest  flexural 
mode  F(0)  and  lowest  symmetric  Lamb  mode  L  (0)  normalized  to  the 
shear  wavenumber  k*  in  the  aluminum.  This  plot  applies  to  both  the 
free  plate,  and  to  the  foam  covered  plate.  Note  that  3  for  the  lowest 
SH  plate  mode  SH(0)  is  equal  to  k'(  . 

The  attenuation  constant  o,  normalized  to  k'  .for  the  leaky  weaves 
is  plotted  in  Figure  41.  Because  of  the  low  values  for  a,  the  waves  can 
propagate  over  several  hundred  wavelegths  before  experiencing  signifi¬ 
cant  attuation.  For  example,  the  SH(0)  mode  at  0.5  MH^  propagates 
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TABLE  V 


Material  Properties  for  Foam  Core  Panels 


Aluminum 

p'(kg/m3) 

2.  70  x  103 

X  (N/ma) 

6.  13  x  1010 

H  (N/m2) 

2.  85  x  1010 

Cg  (m/s) 

3.  25  x  103 

Cp  (m/ s) 

6. 62  x  103 

Foam 

Glue 

48 

5.  8  x  102 

3.  17  x  106 

2.  2  x  109 

L  45  x 106 

5.  6  x  108 

1.  74  x  102 

2 

9.  83  x  10 

3.56  x  102 

2.39  x  103 

-1Q&- 


0.2 


0.6 


1.0 


Figure  40.  Propagation  constants  0  of  the  L(O)  and  F(0)  leaky  wave 

modes  of  a  1/16  in.  aluminum  place  on  which  a  polyurethane 
foam  has  been  formed.  These  curves  are  numerically  the 
same  as  those  found  for  a  free  plate. 
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500  wavelengths,  or  about  3.3  meters  before  attenuating  by  l/e. 

Because  the  foam  can  have  a  high  attenuation  constant  at  those 
frequencies,  we  have  investigated  the  effect  of  loss  on  the  curves  in 
Figures  40  and  41.  No  significant  effect  was  found  up  to  a  loss  tangent 
of  unity.  This  observation  indicates  that  the  rate  at  which  guided  wave 
in  th  aluminum  plate  sheds  energy  into  the  foam  is  not  strongly  depen¬ 
dent  on  the  loss  tangent.  However,  the  foam  attenuation  will  prevent 
the  radiated  energy  from  reaching  the  second  aluminum  plate  in  a  panel. 

B.  Aluminum  -  Glue  -  Foam  Configuation 

The  pretence  of  the  glue  layer  has  two  significant  effects  on  the 
wave  propagation.  First,  even  without  the  foam  it  makes  the  effective 
plate  thickness  greater,  so  that  higher  plate  modes  can  propagate. 

Second  it  acts  as  an  impedance  match  between  the  aluminum  and  foam, 
so  that  the  perturbation  caused  by  the  foam  is  much  greater.  We  first 
consider  these  effects  for  the  SH  modes,  and  subsequently  treat  the 
Lamb  modes. 

1.  SH  Modes 

The  propagation  constant  (3,  normalized  to  k'  is  indicated  in  Figure 

42  for  both  the  aluminum  -  glue  and  aluminum  -  glue  -foam  combinations. 

In  computing  the  curves  of  Figure  42,  the  glue  layer  was  assumed  to  be 

l/64  in.  thick,  which  is  l/ 4  of  the  thickness  of  the  aluminum.  At  low 

frequencies  where  only  the  SH{0)  mode  propagates,  in  the  aluminum  - 

glue  combination  0  is  seen  to  be  nearly  equal  to  k'g,  indicating  that  the 

energy  in  the  wave  is  carried  primarily  in  the  aluminum  with  only  a 

small  fraction  in  the  glue.  As  the  frequency  increases  above  0.6  MH  , 

z 

the  SH(0)  mode  shows  dispersion  and  the  SH(1)  mode  goes  above  cut  off. 
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Figure  42.  Propagation  constants  0  of  the  SH(O)  and  SH(2)  modes 
of  a  1/16  in.  aluminum  plate  covered  with  1/64  in.  of 
glue,  with  and  without  bonding  to  a  semi-infinite 
polyurethane  foam. 


For  the  SH(0)  mode,  0/k'  tends  towards  the  value  3.  3, which  is  the  ratio 

s 

of  the  shear  wavenumber  in  the  glue  to  that  in  the  aluminum.  This  be¬ 
havior  of  0  indicates  that  the  energy  of  the  wave  is  concentrated  in  the 
glue.  The  energy  of  the  SH(1)  mode  for  frequencies  somewhat  above 
cut  off  is  concentrated  in  the  aluminum. 

When  the  foam  is  taken  into  account  the  dispersion  characteristics 

appears  to  shift  to  lower  frequecnies.  Thus  over  the  frequency  range 

computed,  the  SH(Q)  mode  has  0w3.3k*  ,  while  the  cut  off  of  the 

SH(1)  mode  is  below  0.  2  MH^.  For  frequencies  above  cut  off  of  the 

SH(1)  mode,  0  is  intially  close  to  k'  and  then  increases  towards  the 

s 

shear  wavenumber  in  the  glue.  As  before,  this  behavior  suggest  that 
the  energy  of  the  wave  is  carried  primarily  in  the  aluminum  at  lower 
frequency,  and  then  becomes  concentrated  in  the  glue  at  higher  fre¬ 
quencies. 

The  attenuation  constant  a,  normalized  to  k'B,  is  plotted  in  Figure 
43  as  a  function  of  frequency.  For  the  SH(0)  mode,  or  is  larger  than 
that  in  Figure  41  for  the  aluminum  -  foam  conbination  by  two  orders  of 
magnitude.  Thus  l/e  attenuation  length  is  on  the  order  of  five  times  the 
shear  wavelength,  or  5  x  3 . 3  «  17  times  the  wavelength  2ir/0  of  the  guided 
wave.  The  high  attenuation  constant  is  a  relult  of  the  fact  that  the  energy 
of  the  SH(0)  mode  at  these  frequencies  is  concentrated  in  the  glue  layer, 
whose  elastic  properties  are  much  closer  to  the  foam  properties  than 
are  those  of  aluminum.  This  interpretation  is  supported  by  the  SH(1) 
mode  whose  att  enuation  is  lowest  when  the  energy  is  concentrated  in  the 
aluminum  layer,  and  highest  when  the  energy  is  concentrated  in  the  glue. 

From  the  viewpoint  of  NDE,  the  distinct  properties  cited  above  for 
the  aluminum,  aluminum- glue,  and  aluminum- glue -foam  combinations 
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Figure  43.  Attenuation  constants  a  of  the  SH(O)  and  SH(1)  leaky 
wave  modes  of  a  1/16  in.  aluminum  plate  covered  with 
1/64  in.  of  glue  that  bonds  the  plate  to  a  semi¬ 
infinite  polyurethane  foam. 
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facilitate  inpection  of  foam  core  panels  for  regions  of  de- bonding.  Re¬ 
gions  of  the  aluminum  plate  lacking  in  glue,  or  having  glue  but  not 
bonded  to  the  foam  have  very  different  values  of  /3  and  a  ,  and  hence 
can  be  distinguished  by  phase  and  attenuation  measurements.  For  ex¬ 
ample,  consider  a  pair  of  transducers  designed  to  launch  and  receive  SH 
waves  with  j3  =  3.  3k^  at  0,4MHz.  'When  placed  over  a  portion  of  the  plate 
glued  to  the  foam,  the  output  of  the  receiv’ng  transducer  will  have  some 
finite  value  due  to  the  radiation  and  detection  of  the  surface  wave.  Over  a 
portion  of  the  plate  not  glued  to  the  foam,  the  output  will  have  a  lower 
value  since  the  surface  wave  has  a  different  wavelength  and  will  not  be 
strongly  excited  or  detected.  Thus,  unglued  regions  can  be  detected  from 
changes  in  receiver  output. 

2.  Lamb  Modes 

Because  of  coupling  between  P  and  SV  plane  wave  constituents,  the 
dispersion  characteristics  of  the  Lamb  modes  are  more  complex  than 
for  the  SH  modes.  The  propagation  constant  j3,  normalized  to  k'g  is 
plotted  in  Figure  44  for  the  lowest  four  Lamb  modes.  The  labeling  L(0). 
F(0),  F(l)  and  L(l)  is  made  in  accordance  with  their  low  frequency  be¬ 
havior,  since  the  modes  cannot  be  uniquely  separated  on  the  basis  of 
symmetry  in  the  composite  plate.  Curves  for  the  aluminum  -  glue  com¬ 
bination  are  shown  dashed,  while  those  for  aluminum  -  glue  -  foam  are 
solid.  The  dispersion  curve  of  the  L(0)  mode  in  the  presence  of  foam 
could  not  be  traced  past  0.  65  MH  because  of  convergence  difficulties  in 
the  computer  program.  Comparing  the  curves  of  the  L(0)  and  F(0)  modes 
with  those  of  an  aluminum  plate  as  shown  in  Figure  40,  it  is  seen  that  the 
presence  of  the  glue  significantly  influences  the  high  frequency  behavior 
of  the  mode  8. 
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Figure  44.  Propagation  constants  0  of  Lamb  modes  of  a  1/16  in. 

aluminum  plate  covered  with  1/64  in.  of  glue,  with 
and  without  bonding  to  a  semi-infinite  polyurethane 


The  presence  of  the  foam  is  also  seen  to  cause  major  changes  in 
the  propagation  constant,  as  well  as  introduce  attenuation.  In  Figure 
45  we  have  plotted  the  attenuation  constant  a  normalized  to  k'g  for  the 
curves  labeled  F  (0)  L  (0)  and  F(l)  in  Figure  44.  While  the  curves  j3(w  ) 
for  the  three  modes  intersect,  a  for  each  mode  is  different  at  the  inter¬ 
section  frequency,  so  that  all  three  branches  of  the  dispersion  equation 
are  distinct.  The  F(0)  branch  has  high  attenuation  at  low  frequency, 
which  decreases  with  frequency.  Similarly,  the  attenuation  of  the  F(l) 
branch  decreases  dramatically  past  the  cut  off  frequency  of  the  F(l) 
mode  in  the  aluminum  -  glue  combination.  The  low  values  of  a  for  the 

L  (0)  and  F(l)  modes  away  from  0.64  MH  are  consistent  with  the  values 

z 

obtained  in  the  aluminum  -  foam  combination.  The  F(0)  mode  however 
shows  much  higher  a. 

The  curves  of  Figures  44  and  45  indicated  that  the  presence  of  glue 
and  adhesion  to  the  foam  can  be  determined  from  measurement  of  a  and 
fi  for  the  Lamb  modes,  as  discussed  in  the  previous  section.  However, 
a  more  complete  analysis  of  the  mode  properties  and  their  dependence 
on  layer  thickness  is  required  to  devise  effective  test  procedures. 
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Appendix:  Thin  Plate  Approximation  of  the  Plate  Modes 


The  geometry  of  an  infinite  plate  is  illustrated  in  Figure  A- 1 .  The 
total  fields  in  the  plate  can  be  expressed  as  a  superposition  of  the  modal 
fields  of  an  infinite  plate.  In  general,  we  are  concerned  with  modes  propa¬ 
gating  obliquely  with  respect  to  the  x  axis.  These  modes  will  couple  at 
junction  planes  perpendicular  to  x,  at  which  we  must  satisfy  boundary  condi¬ 
tions  involving  the  particle  velocity  v  and  the  surface  traction  x^-T- 

Obliquely  traveling  modes  in  an  infinite  plate  have  x-z  dependence  of 
the  form  exp[- j(+nx  +  k^z)]  where  /k  +  k^  is  equal  to  the  modal  wavenumber. 

In  order  to  satisfy  boundary  conditions  at  junction  planes  perpendicular  to  x, 

all  modes  must  have  the  same  transverse  wavenumber  k  .  Because  of  scattering 

z 

at  the  junction,  plate  modes  having  both  +  x  components  of  propagation  must  be 
accounted  for,  as  indicated  in  the  spatial  dependence  cited  above.  In  addi¬ 
tion,  the  fields  of  the  plate  modes  have  a  variation  with  y,  i.e.,  in  the 
direction  through  the  plate. 


In  the  microwave  network  representation  of  guided  acoustic  waves,  the 
six  quantities  v  and  xq‘T  entering  the  boundary  conditions  are  organized  into 
1x3  matrices  G  and  Q,  which  are 


Here  g(y ,z)  and  £(y,z)  are  the  'hiode  vectors"  describing  the  variation  of 

the  fields  transverse  to  x,  and  the  relative  amplitude  of  the  various 

field  quantities.  As  discussed  above,  the  z  variation  of  g(y,z)  and  £(y,z) 

is  contained  in  the  factor  exp(-jk  z).  The  modal  voltage  V(x)  and  current 

z 

I(x)  satisfy  standard  "electrical"  transmission  line  equations  involving 
the  wavenumber  h  and  modal  impedance  Z,  as  indicated  in  Figure  A-l.  The 
voltage  and  current  are  composed  of  traveling  wave  solutions  having  x  dependence 
exp (+  j nx) . 

r 2  2* 

The  dispersion  relation  giving  the  wavenumber  >/x  +  of  each  mode,  and 

the  transmission  line  representation  of  the  plate  modes  have  been  derived  by 
Markman,  et  al^’^.  In  this  Appendix,  we  approximate  their  results  for  h,  Z 
and  the  mode  vectors  for  the  case  of  a  thin  plate,  whose  thickness  is  small 
compared  to  the  acoustic  wavelength. 

A.  Lamb  Modes  of  an  Infinite  Plate 

Lamb  waves  in  the  plate  are  characterized  by  having  components  of  particle 
displacement  perpendicular  to  the  faces  of  the  plates  and  components  parallel 
to  the  direction  of  propagation. 

1.  Lowest  Symmetric  Mode:  L(0) 

The  symmetric  Lamb  modes  are  those  for  which  the  particle  velocity  in 
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H2-k2\  k 
s  L  s 


kS  ►  4k  A 

-  cot  Kp  -2-  +  -i_£  cot  Ka  |  -  0 


(A-2) 


Here  the  wavenumbers  along  y,  *  and  H  .are  functions  of  k  and  ou  defined 

r  S  L 

by  the  relationships 


k  2  +  h  2  =  k  2  s  u)2  ? — ^ — — 
L  P  P  \  +  2ji 


k T2  +  H  2  =  k  2  a  .2  e. 
L  s  s  w  u 


(A- 3) 


where  X  and  y,  are  the  Lame  constants,  p  is  the  mass  density  and  w  is  the 
radian  frequency. 

The  modal  fields  obtained  by  summing  the  contributions  from  the  P  and 

SV  plane  waves,  which  are  reflected  back  and  forth  between  the  faces  of  the 
plate,  are 


£(y,z) 


—  j  * 

'j  r _  sin  h  y 

/i*P  k  P 


1  k 
1  z 

/fp  kS 


COS  K  y 

p  7 


1  2  _  2  2 
2  A  “  2*  "  ^k  )cos  h  y 

P  Z  '  p  J 

s 


(A-4) 


,  _1_  kL 

J  /TP  \  Sin  %  y 


nP  COB  *D  | 

"S  COS  *S  f 


1  k 
1  Z 

7 =r  cos  h  y 


-j  k  z 
*  * 


'  2  nr  cos  *s  y 

S  L 


-118- 


-119 
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(A- 8) 


$ 


t/2 


can  be  satisfied.  With  this  normalizing  condition,  the  power  carried  along  x 
per  unit  length  along  z,  is  equal  to  the  power  ^  Re  VI*  on  the  modal  trans¬ 
mission  line  for  modes  with  real  or  imaginary  h.  We  impose  this  condition  to 
facilitate  conservation  of  power  test  of  our  results  for  mode  coupling. 

Approximate  solution  to  the  dispersion  equation  (A-2)  can  be  obtained 
for  the  case  of  a  thin  plate  at  low  frequencies  (h  t  «  1  and  k  t  «  11  bv 

8  p  7  1 

approximating  cot  Hp  ~  and  cot  h.  |  as  and  This  approximation  leads 

p  s 

to  a  solution  for  the  wavenumber  k^  of  L(0)  mode,  which  is  the  lowest  symmetric 
Lamb  mode,  and  the  only  symmetric  Lamb  mode  that  propagates  down  to  zero  frequency. 
It  is  found  that 


1  -  v 


'L  J  2  "s  (A- 9) 

where  v  is  Poisson's  ratio. 

With  (A-9)  and  the  approximations  cos  9  «  1  and  sin  0  sw  9  for  0  «  1, 
the  mode  functions  in  (A-4)  and  (A-5)  can  be  approximated  by 
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with  the  characteristic  impedance  chosen  as 


7  =  <“£ 

L  K 

The  normalization  constant  N  obtained  via  (A-8)  is  given  by 

Li 
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2 .  Antisymmetric  Modes:  F(0)  and  F(l) 

In  the  antisymmetric  modes  v  is  an  odd  function  of  y  whereas  v  is  an 

x  y 

even  function,  with  the  result  that  the  plate  bends  or  flexes  during  its 
motion.  These  modes  are  therefore  referred  to  as  flexural  modes.  The  dis¬ 


persion  relation  giving  the  wavenumber  k 

F 

function  of  frequency  is 
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where  the  wavenumbers  Rg  and  Kp  are  defined  by  the  relationships 
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The  modal  fields  are  obtained  again  by  summing  the  fields  contributed 
from  P  and  SV  plane  waves  that  are  multiply  reflected  between  faces  of  the  plate. 
The  mode  vectors  are  given  by 
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and  n  are  given  in  (A-6)  with  k  replaced  by  k  . 
p  L  c 

characteristic  impedance  of  the  modal  transmission  line  is  chosen  as 
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where  k  =  / k  -  k  is  the  wavenumber  along  x.  The  signs  of  q  and  ?  are 
*r  F  z  F 

chosen  such  that  (A-8)  is  satisfied. 

There  are  two  significant  flexural  modes  that  should  be  taken  for  the 
wave  propagation  in  the  thin  plate  at  low  frequencies.  One  of  these,  the 
lowest  flexural  mode,  F(0),  has  a  real  wavenumber  even  for  frequency  approach¬ 
ing  zero,  while  the  other,  the  first  higher  flexural  mode,  F(l),  has  an 
imaginary  wavenumber.  The  F(l)  mode  cannot  be  neglected  even  though  it  has 
an  imaginary  wavenumber  since  the  imaginary  wavenumber  vanishes  as  *-he  fre¬ 


quency  goes  to  zero,  and  since  its  transverse  dependence  is  similar  to  that 
of  the  F(0)  mode.  As  a  result,  the  presence  of  F(l)  strongly  influences  the 
amplitude  of  the  F(0)  mode,  and  other  modes,  excited  at  a  junction.  In  addi¬ 
tion  the  F(l)  couples  neighboring  junctions  at  low  frequencies. 


The  wavenumbers  of  the  lowest  flexural  mode,  and  the  first  higher 


flexural  mode,  k  , can  be  obtained  from  the  dispersion  equation  (A- 14)  by 

F1  t  V  V  V 

approximating  tan  *  —  and  tan  — —  as  — jj—  and  — for  the  case  of  a  thin  plate 


at  low  frequencies.  The  modal  wavenumbers  are  found  to  be 


(A- 19) 

The  mode  functions  given  in  (A- 16)  and  (A- 17)  may  be  approximated  by 

e2 

substituting  sin  9»9  and  cos  9  <=«  1 - which  is  appropriate  for  the  case 

of  a  thin  plate.  The  thin  plate  approximations  for  the  F(0)  mode  gives 
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and  the  characteristic  impedance  is  given  by 
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The  normalization  constant  can  be  obtained  via  (A-8)  as 
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The  thin  plate  approximations  for  the  F(l)  mode  gives 
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where  r  is  given  in  (A-22) .  For  the  F(l)  mode 
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The  -  sign  in  the  characteristic  impedance  and  the  mode  functions  q  has 
been  chosen  in  order  to  satisfy  normalizing  condition  (A-8). 

B.  SH  Modes  of  an  Infinite  Plate 

The  SH  modes  are  orthogonally  polarized  with  respect  to  the  Lamb  modes 
in  that  the  particle  displacement  of  an  SH  mode  is  parallel  to  the  plate 
faces  but  perpendicular  to  the  propagation  direction.  The  lowest  mode, 
SH(0),  is  of  even  symmetry  and  the  only  SH  mode  that  propagates  down  to  zero 
frequency.  The  dispersion  relation  for  symmetric  SH  modes,  which  yield  SH 


wavenumbers ,  k  ■ 
SH 

simple  expression 


2 


+  k^  as  a  function  of  frequency,  is  given  by  the 


kgH  «/kg2  -  (2nrr/t)2  n  -  0,  1,  2, 

Here,  k  *  w  /®-  is  the  shear  wavenumber, 
s  VP- 


(A-29) 


The  normalized  mode  functions  of  even  symmetry  are  given  by 
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where 
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The  characteristic  impedance  is  chosen  to  be 

Z  =  M  (A- 34) 

S  0) 

The  dispersion  relation  and  the  mode  functions  for  the  lowest  SH  mode, 
SH(0),  can  be  obtained  by  simply  setting  n  =  0  in  (A-29)  through  (A-33). 

One  obtains 

kSH  '  ",  <A‘35) 
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The  modal  quantities  used  In  the  teat,  and  Hated  In  equations  (7)-(U) 
ate  obtained  by  setting  t,  *  0  *«  cotte.ponding  expression  Hated  in  this 

Appendix. 
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